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ABSTRACT 

Descriptions of several digital computer programs for electronic 
circuit analyses are given. A listing of a special-purpose (Class M C n 
amplifier design) program is given. Comparisons of several general 
purpose programs are made. An extensive description of the CADLIB 
is included. The discussion of the CADLIB points out numerous 
modifications to the general purpose programs. These modifications 
include such things as the addition of printer-plots, more accurate 
algorithms, etc. Also many of the problems encountered in adapting 
the various programs to the NPGS machine (IBM 360 model 67) are 
enumerated. Appendicies are included to provide a reference for coding 
the members of the CADLIB for which no instruction manual was avail- .. 
able. An extensive bibliography on applied computer aided design and 
the job control language details of how the CADLIB was created are 
included. In addition two extensive examples using the CADLIB are 



included . 
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CHAPTER I 



INTRODUCTION 

In circuit design a most time-consuming, error-prone, and 
frustrating chore is the repetitious solving of network equations or 
the repetitious measurements taken from a breadboarded circuit as 
some parameter(s) is varied. However, the digital computer thrives 
on simple-minded repetitious work while attaining perfection in the 
writing of equations and near perfection in their solutions. The 
accuracy of the solutions may be even further enhanced by the use of 
double precision calculations which, under the present trend of com- 
puter design, is very nearly as fast as single-precision calculation, 
e.g. , the ADD operation in the IBM 360 requires 2.43 microseconds 
in single-precision and 2.45 microseconds in double-precision. 

Another point in regards to accuracy is that no solution is any 
more accurate than its mathematical model irrespective of the number 
of significant digits carried in the computation process. Modelling 
of resistors, coils, and capacitors with the necessary fringing or 
stray effects presents few problems. However, an accurate model 
for an MOS-FET or integrated circuit is prohibitively complicated 
for the slide- rule or calculator approach. But by making use of the 
speed of the computer virtually any model can be effectively used. 

In the following discussions a basic knowledge of the Fortran 
programming language is assumed. There are also numerous 
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references to OS/ 360 job control language but the understanding of 
these references is not necessary to make use of this report. In 
particular the discussions on how the various programs were placed 
on a disk is intended primarily for guidance to one who is interested 
in modifying or adding to these programs. 
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CHAPTER II 



COMPUTER AIDED DESIGN PROGRAMS 
INTRODUCTION 

There are several hundred digital computer programs presently 
in use for circuit design. These programs run the gamut from simple 
specialized programs for a particular circuit configuration to extremely 
complicated packages that have widely varying capabilities such as: 
sensitivity coefficients, worst case, DC, AC, and transient analysis, 
plotted outputs, and recently on-line graphics (light pen and CRT). 
Appendix A contains an example of a specialized program to "optimize" 
a Class "C" amplifier feeding a resistive load. As is characteristic 
of these specialized programs the Class "C" program is of limited 
value and the input is very restrictive both in content and form. Several 
worthwhile specialized programs in various Fortran languages have 
been listed and documented in JPL ! s Technical Report 32-979, 
references 18 and 19. Examples of the general purpose type of 
program comprize the remainder of this report. 

GENERAL DISCUSSION 

As pointed out above the capabilities of the available programs 
vary so widely in applicability and scope that any comparison must be 
limited to the few large general-purpose programs. Table 1 lists the 
major characteristics of some of the major programs. The input to 
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these programs is usually engineer- oriented and free-format , i.e. , 
not Fortran. Typically the user is required to: 

(1) . Number the nodes sequentially from zero or one, 

(2) . Either number the branches sequentially or give them a 

unique name, and 

(3) . Code the branches on cards using commas as delimiters. 

The limitations of these programs must also be considered. 

All state- variable based programs suffer from two common limitations. 
These are: 

(1) . Numerical errors including round-off errors and 

truncation errors , and 

(2) . The time constant problem generated by the network 

being analyzed. 

The first limitation is caused by the limited number of digits in a com- 
puter word and can be minimized by the use of better algorithms and 
double-precision arithmetic. The second limitation is the result of 
solving continuous integral equations by quantized numerical methods. 
To obtain accuracy of integration a small step size is needed for even 

"'’By engineer- oriented we mean that a branch is described by its 
node numbers and a positive direction of current flow. 

2 

By free format we mean that specific columns on the card are 
meaningless and delimiters such as commas are used to separate the 
data. 
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ECAP 


NET-1 


CIRCUS 


SCEPTRE 


PREDICT 


NASAP 


Init. Cond. 


X 


X 


X 


X 


X 


X 


Sens. Coeff. 


X 






X 


X 


X 


DC Analysis 


X 


X 










AC Analysis 


X 












Trans. Anal. 


X 


X 


X 


X 


X 


X 


Max Element 


200 


600 


400 


300 


300 




Max Nodes 


50 


100 


200 


300 


100 




Stored Model 




X 


X 


X 






Volt. Sou. 


X 1 ' 2 


x 1 ’ 2 


x 1 ' 2 


x 1,2,3 


X 1 - 2 




Curr. Sou. 


X 1 - 2 


x 1 ’ 2 


X 1 - 2 


x 1 ’ 2 ’ 3 


x 1 ’ 2 




Language 


FTN 


FAP 


FTN 


FTN 


FTN/ F AP 




Deck Size 


5700 


30,000 


10,000 


15,000 


10,000 




F ormulation 


Nodal 


St Var 


St Var 


St Var 


St Var 


Dich 


Integration 


Imp . 
Numer . 


Pre. 
Corr . 


Exp . 


Exp , Trap 
& R-K 4 


Trap 





Notes: Can be constant 

2 

Can be function of time, table, or equation 
3 

Can be a function of a subroutine 



Table 1. A comparison of the Major Programs. 
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the smallest time constant of the circuit. But to obtain the complete 
transient or steady- state solution the integration must be carried out 
to about five of the largest time constants of the circuit. Thus we once 
again hit that perennial "trade-off" problem and about the only solution 
at the present is to eliminate as many of the small time constants as 
practical and/or use a variable step-size integration technique to 
minimize the error. 

Another method of circumventing the first limitation is to use 
a topological approach. However, with growing circuit size computer 
run time rapidly becomes excessive due primarily to the requirement 
that all possible trees must be found. This tree-finding problem will 
be treated more thoroughly in a later section dealing with the CALAHAN 
program. At the moment topological methods are limited to networks 
of 25-30 branches and 15-20 nodes. 

ECAP (Electronic' Circuit Analysis Program) 

ECAP [16] was developed through a joint effort of IBM and 
the Norden division of United Aircraft Corporation and released in 
June 1965 for the IBM 1620 computer. It has since been translated for 
use on the IBM 7040/44, 7090/94, and the System 360 series. ECAP 
with modifications is probably the most widely used general purpose 
CAD (Computer Aided Design) program today. It is one of very few 
(if not the only) general purpose programs to offer both DC and AC 
analysis with an automatic parameter modification feature in both its 
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DC and AC analyses. ECAP will provide, as output, the sensitivity 
coefficients, worst case, and standard deviations all of which are 
noticeably lacking in many other general purpose programs. Another 
positive feature is the switch (not available on most other programs) 
which provides a very flexible piece-wise linear modelling capability 
to approximate non-linear elements. 

A major limitation of ECAP is its lack of a non-linear capability 
except for its extremely good piece- wise approximation capability. 

The major problems in using ECAP are its lack of a model library for 
frequently used components and the problem of data processing, e.g. , 
a single node voltage is unobtainable but must be separated from a 
tabulated list of all node voltages. 

NET-1 

NET-1 £lo3 was developed at Los Alamos Scientific Labora- 
tory of the University of California at Los Alamos, New Mexico, under 
the auspices of the United States Atomic Energy Commission. NET-1, 
like ECAP, is a general purpose program but suffers most noticeably 
from the lack of an AC capability and the unfortunate (?) choice of 
FAP as a primary programming language; thereby, limiting its users 
to the IBM 7090/94 and Maniac II computers as conversions to other 
machines would be a major undertaking. However, NET-1 does 
have a built-in non-linear model for transistors and diodes but these 
models require some 36 parameters for the transistors and 13 for 
the diode and the model must be entered before it can be used. These 
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parameters are not readily available from manufacturers and must be 
measured from several devices of each type to obtain typical values. 
NET-1 became available in October 1962 and has since been rather 
widely distributed with its users swapping semiconductor character- 
istics rather freely. 

CIRCUS 

The CIRCUS program \j.Z & 13^ was developed at the Boeing 
Company, Seattle, Washington, primarily for use in studying the 
effects of nuclear radiation on electronic circuits, but its applicability 
to transient analysis of all types led to its use as a general purpose 
program. CIRCUS has the modeling capability (and limitations) of 
NET-1 but was written in Fortran IV (98%) and Assembly (2%) languages, 
and the conversion to the System 360 has been made by Lincoln Labora- 
tories at MIT. CIRCUS is relatively fast and accurate but suffers 
from the same fixed models as NET-1. This can lead to unnecessarily 
long run times due to the small internal time constants of some charge- 
control models. 

SCEPTRE 

Like CIRCUS a desire for a transient analysis of the effects 
of nuclear radiation led to the development of SCEPTRE [h] . it 
was developed by IBM for Kirkland AFB. One of the contract stipula- 
tions was "100% Fortran" and this led to a very machine dependent 
program completely defeating the purpose of the Fortran restriction. 
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SCEPTRE has two unique features not found in other programs. 

The first is the ability to scan the loop matrix of the circuit, 
select an efficient partition, generate a Fortran program describing 
the circuit, and then execute the generated program. If desired, this 
Fortran program can be saved and executed at a later time or in steps 
as might be desired when the run time becomes long. Parameter 
modifications make changes only in this generated program so that the 
load time (time spent generating the Fortran program) is encountered 
only once for each circuit. This load time originally required 4-5 
minutes of run time; however, this has been lowered to about 2 minutes 
[Zl] through the use of more efficient overlays. 

The second unique feature of SCEPTRE is its tremendous 
modelling flexability. SCEPTRE will accept and store any circuit in 
its model library. The connections and their order is specified at 
the time of entering the model into the library and a simple call will 
connect the model anywhere it is desired. The value of the model 
library for ’’standard 11 or often used circuits is readily apparent. Due 
to its size and its modelling capability SCEPTRE is probably potentially 
the most powerful general-purpose CAD program available today. 

PREDICT 

PREDICT is the forerunner of SCEPTRE but as originally for- 
mulated it has no stored model capability and had a rather simple- 
minded trapozoidal integration algorithm. However, an IBM 360 



17 



version has been written and some of the features of SCEPTRE were 
added to PREDICT to form PREDICT- 360. 

NASAP (Network Analysis for System Application Program) 

NASAP £4 & 7^ is mentioned primarily because it has had a 

maximization of advance publicity in spite of the fact that it is not 

operational on any machine at the moment. But NASAP is planned to 

have virtually all of the features of most of the available general-purpose 

programs available today. NASAP requires no lengthy matrix inversions, 

thereby eliminating the primary source of round-off and truncation 

errors. NASAP utilizes a dichotomous approach with a topological 

analysis of many small circuits and preliminary testing indicates that 

5 

the computer run time is reduced by factors as large as 10 from the 
state variable based programs. 
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CHAPTER III 



CADLIB 

The CADLIB is a computer library of computer programs 
designed to aid the designer in the analysis of various circuits. The 
CADLIB is presently maintained in executable form on a direct access 
disk named J0BOO3 in the computer center of the Naval Postgraduate 
School. Support for the CADLIB is furnished in the form of a sequential 
data set named MODEL (also on J<Z>B003) for use as a model library for 
one of the members of CADLIB. The various members of CADLIB are 
executable members of a partitioned data set and they consist of the 



following: 




ECAP 


slightly modified by J. G. Sketoe to furnish 
printer plots £l73 . 


CIRCUS 


OS/ 360 version distributed by MIT except 
that the default model unit 1 1 capability 
was inserted. 


CALAHAN 


slightly modified by J . G. Sketoe to furnish 
automatic printer plots both of freq and 
transient responses. 


POTTLE 


slightly modified by J. G. Sketoe to furnish 

automatic printer-plots of the freq response. 
Also the matrix inversion has been changed 
to double-precision arithmetic. 


POTTOL 


the POTTLE program with a free-field format 
capability as designed by J. G. Sketoe 
(See Appendix D) . 


MODEL 


the CIRCUS modeling program as described in 
reference 13. Operates on data set MODEL. 
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The CADLIB data set was constructed using a SYSLM0D of 
the following form: 

/ / SYSLM0D DD DSN AME= CADLIB , UNIT=2311 , DISP=(NEW,KEEP) , X* 

// V<Z>LUME=SER=J0BOO3, 

// SPACE=(CYL, (25 ,2,2)), 

// LABEL=EXPDT=99350 

And the individual members were added using 

/ / SYSLM0D DD DSN AME= CADLIB , UNIT=2311 , DISP=0LD, X 

V</>LUME=SER=J(/)B003 

Execution of the various members requires different sequences 
of control cards. These sequences will be given below as each member 
is discussed. 

ECAP 

The first program obtained for the CADLIB was IBM's Elec- 
tronic Circuit Analysis Program. The 360 version requires approx- 
imately 290K bytes of core when used without overlays as in CADLIB. 
ECAP was obtained from Share, a facility of IBM's, in the E-level 
compiler edition. As the OS/ 360 version was ordered only the 360 
manual [Z0] was included. However, this manual only pointed out the 
differences in ECAP-1620 and ECAP- 360. A manual was then ordered 
for the 1620 version p6] and two months later it arrived. In the mean- 
time by using ECAP problems published in various journals it was 

This X as in all of the following control cards must be in column 72. 
Also the first character of the following card (except for the double 
slashes) must be in column 16. 
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discovered that the received version simply would not work properly. 
The major problems and their solutions were: 

(1) . The release 9 G-level "in-house* 1 compiler failed to con- 

struct the proper instructions for overlay. At that time 
the computer center policy was to run two-256K byte 
machines rather than one 512K byte machine so that 
various hardware testing could be done concurrently with 
computer service. Use of the older but slower H-level 
compiler cured this trouble. (This fault of the G- compiler 
has been cured in later releases). 

(2) . The next problem discovered was that the language 

processor portion of ECAP did not function correctly 
when compiled on either the G- or H-level compiler. The 
cause of this trouble lay at the fundamental differences 
in the E-level Fortran compiler and the G/H level com- 
pilers. As the language processor converts an engineer- 
oriented circuit description into the required computer 
description, literal data (alphanumeric characters) is 
manipulated vice purely binary numbers (floating point 
numbers are a special coding of binary). This manipula- 
tion is done by using a unique combination of bits to 
represent a particular character. When an arithmetic IF 
statement is used to decode this data, the data is treated 
as numbers and not as literal characters. For example, 
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the Fortran statement 



IF (M-N) 10, 15, 10 

compiles into an arithmetic subtraction of N from M, 
examine the resultant, then branch to statement 10, 15, 
or 10 according to whether the resultant is negative, zero, 
or positive. The contents of a computer word or cell is 
considered positive or negative according to whether the 
first bit is a zero or one, respectfully. In the IBM 7090/94 
machines the literal data codes were mixed as positive 
and negative numbers with the special characters (, . etc.) 
intermingled within. This mixture was observed in the 
E-level machines while in later machines this coding was 
changed such that alphanumerics are negative (and 
sequential) and special characters are positive. The 

solution to this problem was quite obvious once the source 

1 1 ■ 

of the error was known. All literal data arithmetic IF's 
were replaced with logical IF*s which do their comparisons 
on a bit-for-bit basis rather than an arithmetic answer. 

(3), Solving of these two problems finally got our copy of ECAP 
at least partially working. The next problem encountered 
was that the automatic parameter variation feature did 
not function properly. A cry for help to IBM was answered 
with a list of errors and omissions in the originally dis- 
tributed version of ECAP. The appropriate entries cured 
this problem . 
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(4). Next an addressing error was encountered when trying 
to obtain a frequency response solution with an overlaid 
program. As the various control sections in overlays 
complicate debugging, an attempt was made to find the 
error while using all 512K bytes of core without overlays. 
The addressing error disappeared. Reinsertion of the 
overlay brought the bug back. The only conclusion was 
that the bug was in the machine constructed overlay. 

Next run time comparisons were made on overlay vs. 
non-overlay versions. The results were that overlaid 
programs required 2-4 times as much run time as non- 
overlay. By this time all 512K was being used for almost 
all runs and the decision to use ECAP without overlay in 
the CADLIB was made. The bug was brought to the attention 
of the computer center personnel and left there. 

Presently all known bugs have been worked out of the NPGS 
version of ECAP in addition to the addition of a printer plotting addition 
as explained in reference 17. As User’s Manuals are in short supply 
an abbreviated manual has been written and is published as reference 17. 

The required control cards to execute the CADLIB version of 
ECAP are shown in figure 1 . 
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//ECAP J0B 


(EACb , 0 1FT) , ' NAME ' , MSGLEVEL= 1 




//JOB LIB DD DSNAME=CADLIB,UNIT=2311,DISP=0LD, 


X 


// 


VOLUME=SER=J0BOO3 




// EXEC PGM ECAP 




//FTO6F001 


DD SY SOUT= A 




/ / FTO5F001 


DD * 




Batched 


ECAP coded input decks 




/* 







Figure 1. Control Cards to Execute ECAP. 



Basically ECAP is formulated to fill the nodal conductance 
matrix as the input branches (B- cards) are read in. The various 
command cards then set the appropriate flags to indicate which 
algorithms are executed to obtain the desired outputs. This conductance 
matrix is available as an output request under the MISCELLANEOUS 
request. 

One additional note worthy of comment about ECAP is that it 
only interprets the first two characters of any word, e.g. , PR is the 
same as PRINT, V0 is the same as V0LTAGES, BE is the same as 
BETA, etc. 

CIRCUS 

A. copy of the CIRCUS- 360 program was obtained from Lincoln 
Laboratory, MIT, Lexington, Massachusetts. As with ECAP the 
received version would not work although MIT had exactly the same 
type machine as NPGS, i . e . , a model 65/ 67 IBM 360 . The MIT des- 
cription of the distributed CIRCUS tape is given in Appendix B. The 
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various problems encountered with our copy were: 

(1) . The CIRCUS tape, as received, had an H-level compiler 

output version on file # 1 (Appendix B) which should have 

loaded directly from the tape. However, a tape error 
❖ 

had occurred and the tape had to be recompiled. 

(2) . At MIT there exists three resident routines that access 

the day, year, hour, minute, and second. These 
routines were used by CIRCUS- 360. As these routines 
did not exist at NPGS the initial solution was to provide 
dummy routines to get CIRCUS operating. While this 
solution worked it eliminated the variable time limit of 
CIRCUS as explained in the User ! s Manual [ 12 ] . The 

author then proceeded to learn enough assembly language 
to write a set of routines to accomplish this. These 
routines are now added to our version of CIRCUS as well 
as entered in the resident library of our computer center. 

(3) . While translating CIRCUS MIT evidently inserted a de- 

bugging print-out at each integration step. As this led 
to excessive and useless print-out of the integration step 
time and size, both in time and paper this routine was 
slightly modified to eliminate this waste. 

A bEND card had become BEND where "b" implies a 12-2-9 punch. 

A H B M character implies a 12-2 punch. 
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(4) . No overlay table was included in CIRCUS- 360 as it will 

fit into a 256K byte machine if the load point is below 
B000 (See Appendix B) • However, the load point at NPGS 
is B820jk with just the M in-house M G-compiler thereby 
necessitating an overlay. This overlay was constructed 
and is shown in figure 2. This overlay was used but was 
soon dropped for the following reasons: 

(a) . Substantial decrease in run time. 

(b) . Prior experience with ECAP indicated 

errors might occur. 

(c) . Using the 256K byte machine only 158K 

was available with HASP. CIRCUS even 
with this overlay required 153K leaving 
almost no room for increased compiler 
size (which came later). 

(d) . By going to 512K we were able to increase 

our maximum circuit size to circuits of 
over 70 equations (Appendix B) . 

(5) . The MIT description lists two known bugs (Appendix B) 

in their program. The second one of these failed to show 
up in our recompiled version of CIRCUS. This led to a 
successful attempt to correct the tape error mentioned 
earlier so that this bug could be checked out. As indicated 
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Root 

Segment 


MAIN 

CLOCK 

ENDJOB 

EQUATI 

KLOCK 

KLOCKI 

MATOUT 

MONKEY 

OUTMAT 

TIGER 

COMPR 


PAC 

UNPAC 

EXTRCT 

GET 

PUT 

BCDFPT 

SETLMT 

BCDINT 

TIMHMS 

TIMHR 

TIMDAT 


Overlay 










Alpha 


CHANGE 




MAIN 2 






CREATE 




DOUBLE 






DICT 




MPLY2 






ELFIND 




PVALUE 






FETCH 




SOURCE 






GENRAT 




SVALUE 






INCIDN 




TAB LIZ 






JFIND 








Overlay 


KFIND 










LINKI 


BEAR 


CHARE T 


FINDER 


Beta 


LINK2 


ELPHNT 


EXPINT 


LINK 6 




MAPPER 


HORSE 


GAMT 


PLOTER 




PLSTOR 


INDUCT 


INTRPT 


CLEAR 




REGIT 


LINK 4 


IPPT 


PLACEI 




SETUP 


TMPLY 


LINK 5 


PUTPL 




SHIFT 


TRIPLE 


LISTVL 


FINDY 




SMODEL 




NODEV 


SCALEY 




STORY 




OUTPUT 


FINDX 








PHOTO 


SCALEX 








PLOT 


SPTPLT 








PRINT 


PLTOUT 








PTIME 










SPCAL 










SPCATD 










VMPLY 










XXMIN 







Figure 2. CIRCUS- 360 Overlay Requiring 
approximately 153K bytes. 
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in Appendix B the MIT version failed to recognize a 
card of the form 

USE MODEL LIBRARY, n 

where n is integer denoting the logical unit of the model 
library. The NPGS version did recognize the card 
correctly. Naturally the NPGS version was placed in 
CADLIB . 

(6). An error of omission by MIT but not mentioned by them 
in their write-up was that the default to unit 1 1 (page 2,5 of 
User 1 s Guide Q2] ) was omitted. This has been corrected 
in the NPGS version. 

Presently CIRCUS- 360 is in CADLIB and is working properly. 
As CIRCUS uses several different logical units a cataloged procedure 
has been established for CIRCUS in order to simplify the job control 

i 

card requirements. The standard execution of CIRCUS using logical 
unit 11 for the model library, when needed, is shown in figure 3. 



//CIRCUS J0B (EAC6,01FT),' NAME 1 , MSGLEVEL= 1 
//JOB LIB DD DSNAME= CADLIB , UNIT=231 1 , DISP=0LD, X 
// ' VOLUME=SER=J0BOO3 

// EXEC CIRCUS 
/ / GO . SY SIN DD * 

Batched CIRCUS coded input decks 

/* 



Figure 3. Standard CIRCUS Job Control Language. 
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If logical assignments other than or in addition to the default 
assignments are desired (such as use of "SAVE TAPE"), these units 
must be defined in a manner similar to that shown in figure 4 where 
unit 15 is defined for use as a "SAVE TAPE 1 or as a model library tape. 

A listing of the actual cataloged procedure to execute CIRCUS 
is given in figure 5 for reference althouth it should not be needed 
except to perhaps extend or replace it. 



//CIRCUS J0B (EAC6.01FT),' NAME ' , MSGLEVEL= 1 

//J0BLIB DD DSNAME= CADLIB , UNIT=23 1 1 , DISP=(/)LD 3 X 

// V<Z>LUME=SER=J(Z>B003 

// EXEC CIRCUS 

//G0.FT15FOO1 DD DSNAME=SAVE, UNIT=2400 , DISP=(NEW, KEEP) , X 
// V0LUME=SER=XXX,LABEL=YYY, X 

// DCB=(BLKSIZE= 1204, RECFM=V) 

//G0.SYSIN DD * 

Batched CIRCUS coded input decks 

/* 



Figure 4. Defining Logical Unit 15 for use by CIRCUS 



MEMBER NAME 


CIRCUS 




//GO 


EXEC PGM= CIRCUS 




/ / FTO1F001 


DD DSNAME=SYS1.SYSUT1,DISP=0LD, 


X 


II 


DCB=(BLKSIZE= 1204, RECFM= V 




// FTO3F001 


DD DSNAME=SYS1 . SYSUT2 , DISP=0LD, 


X 


II 


DCB=(BLKSIZE 1 204 , RECFM= V) 




II FT12F001 


DD DSN AME= SY SI . SYSU T 3 , DISP= 0 LD , 


X 


II 


DCB=(BLKSIZE= 1204, RECFM= V) 




/ / FT11F001 


DD DSNAME=M0DEL,UNIT=2311, DISP=0LD, 


X 


II 


VOLUME=SER=JOB003 




//FT06F001 DD SYS0UT=A 




/ / FT05F001 DD DDNAME=SYSIN 





Figure 5. Cataloged Procedure "CIRCUS" 
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The formulation of CIRCUS is based in the state varible approach 
as advanced by Bryant [*22*] and extended by Wilson and Massena [23-] . 
A detailed, although not rigorous, explanation of the mathematical 
theory of CIRCUS is given in the Program Manual fl3^ • 

MODEL 

The MODEL member (alias UPDATE) of CADLIB is a program 
to list, modify, generate, etc. the model library as described in 
Chapter 4 of the CIRCUS Program Manual £l3^J • The CIRCUS tape 

as received from MIT had a compiled version of the MODEL program 
as file #2 (see Appendix B) . This file was used as input to the linkage 
editor to form member MODEL which operates on the actual CIRCUS 
model library. The control cards required for execution of the member 
MODEL is shown in figure 6. The model library is described by the 
FT11F001 card in figure 6 which, as shown, describes a date-protected 
sequential data set named MODEL which resides on the same JOB003 
disk as the CADLIB. If one wishes to use the program or member 
MODEL to operate on some other data set (i.e. , magnetic tape) only 
this FT11F001 must be changed. 
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//MODEL 


JOB (EAC6.01FT), ' NAME 1 , MSGLEVEL= 1 


//J0BLIB DD DSN AME= CADLIB , UNIT=231 1 , DISP=0LD, X 


// 


VOLUME=SER=JOB00 3 


// EXEC 


PGM= MODEL or UPDATE 


//FT11F001 


DD DSNAME= MODEL, UNIT=2311,DISP=0LD, X 


II 


V0LUME=SER= J0BOO 3 


//FT12F001 


DD DSNAME=SYS1 . SYSUT1 , DISP=0LD 


//FT06F001 


DD SYS0UT= A 


//FT05F001 


DD * 


Model input as per Chapter 4, CIRCUS 


/* 


Program Manual [l3J 



Figure 6. Use of the Model Program. 



POTTLE 

The copy of C. Pottle* s Linear Network Analysis Program as 
received required only changes in the logical unit assignments to make 
it compatible with the NPGS machine. To make the use of this program 
feasible from the CADLIB all outputs were assigned to unit # 6 (printer) 
and the input was assigned to unit #5 (card reader). However, flexibility 
was added to the program by making the unit designations variables 
in the individual read and write statements and then assigning these 
variables in a central location in the main program. These variables 
are then transferred between various subroutines via labeled common. 
This provides the user with the facility of writing out some data 
separately from other data. The variable names and their assignments 
are shown in figure 7. 
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Name 


Unit 


Use 


INTAPE 


5 


User’s input data 


IOUTP 


6 


Copy of data set & errors 


STATP 


6 


State equations 


TRANTP 


6 


Transfer functions, poles, & zeros 


FREQTP 


6 


Frequency response 



Figure 7. Default POTTLE I/O Assignments. 



The control cards required to execute the POTTLE program 
from the CADLIB is shown in figure 8. 



/ / P0TTLE J0B (EAC6 , 0 1FT) , !l NAME 1 , MSGLEVEL= 1 




// J0BLIB DD 


DSNAME= CADLIB , UNIT=2 31 1 ,DISP=0LD, 


X 


// ■ 


V<J LUME=SER= J0BOO3 




// EXEC 


PGM=P0TTLE 




//FTO6F00I 


DD SYS0UT= A 


, 


/ / FT05F001 


DD * 




POTTLE coded input (See Appendix C) 

/* 





Figure 8. Control Cards to Execute POTTLE. 



The procedure ^143 implemented by Mr. Pottle is based on a 
recursive algorithm developed by Derisoglu [.6^ which extends Bryant*s 
algorithm to nonlinear active elements. The state equations are formula- 
ted; then the frequency response is determined via the Souriau-Frame 
algorithm to obtain the matrix (sI-A^. The time responses are 
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obtained by factoring the characteristic polynomial and performing 
partial fraction expansions of each element of the transfer function 
matrix or by direct integration of the set of first-order equations. 

The instructions for coding the POTTLE program are included 
as Appendix C. 

Changes to the POTTLE program were limited to (1) providing 
extended accuracy through the use of a double-precision vice the 
original single-precision matrix inversion routine and (2), the addition 
of an automatic plot of the time response and the linear frequency 
response . 

POTTOL 

The frustrations encountered in using the POTTLE program 
pointed out the desirability of a translator so that a free-format input 
could be used. The author has designed such a translator and combined 
it with the POTTLE program described above to form the POTTOL 
program. While developing the Translator several other convenience 
features were included. 

While the use of the POTTOL input format is relatively simple, 
a complete description of the Translator is not so simple. The complete 
description of the Translator input coding is given in Appendix D. The 
control cards to implement POTTOL are shown in figure 9 and are 
very similar to those used for ECAP and POTTLE. 
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// P0TT0L J0B 


(EAC6.01FT),' NAME ' , MSGLEVEL= 1 




//J0BLIB DD DSNAME=CADLIB , UNIT=231 1 , DISP=(#LD , 


X 


II 


V<Z>LUME=SER=J0BOO3 




II EXEC 


PGM=P0TT0L 




//FT06F001 DD 


SYS($UT= A 




/ / F T05F00 1 DD 


* 




Batched POTTOL coded input (See Appendix D) 




f «T“ 







Figure 9. Control Cards to Execute POTTOL. 



CALAHAN 

Probably the most accurate program in the CADLIB is the 
CALAHAN [3*3 which was developed by D. A. Calahan while at the 
University of Illinois. As received, the program was coded in IBM 
7090 Fortran II. The problems encountered with the program were 
generated primarily by a complete lack of knowledge of the intricacies 
of Fortran II. But the differences were not as great as feared. The 
names of the Fortran library routines were slightly different (e.g. , 

COSF vs. COS) and the passing of arguments through subroutines could 
be handled differently causing the computer to get lost when passed 
in the same manner in Fortran IV. The major problem encountered 
in implementing and testing the program was the intepr etation of the 
instruction manual. Appendix E is a set of instructions which (hopefully) 
will clearify the coding of the input to the CALAHAN program. 

The CALAHAN program is a topologically based program which 
is the primary source of its accuracy. This approach evaluates the 
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determinants and cofactors of the node admittance matrix by finding 
the related n-tree products of the corresponding linear graph. Basically 
the matrix solution for the node voltages is transformed by the Binet- 
Cauchy theorem | 1 5 ] and the Maxwell 1 s formula to a problem of 
determining all tree and associated n-tree products. This tree-finding 
problem quickly replaces the integration time constant problem of the 
state- variable program as a computer run time limitation. Regardless 
of the tree-finding scheme used the possible combinations of branches 
increase factorially (about the most rapidly increasing function known) 
and the number of trees increase exponentially (about the second most 
rapidly increasing function known). For example, the circuit of figure 
10 using the hybrid-pi transistor model results in a circuit containing 
eleven nodes and twenty-eight branches which forms 456, 172 trees! 

While all topologically based programs have this tree-finding limitation, 
they do not suffer from the round-off error so prevalent in the state- 
variable formulations. 




Figure 10. Wide-band R-C Coupled Amplifier. 
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The execution of the CALAHAN program from the CADLIB is 
shown in figure 11 and again the similarity to the execution of the 
ECAP and POTTLE should be noted. 



/ / CALAHAN 


JOB (EAC6.01FT, ' NAME ',MSGLEVEL=1 




/ / J0BLIB DD 


DSNAME= CADLIB, UNIT=231 1 ,DISP=(Z>LD, 


X 


// 


VOLUME- SER=J0BOO 3 




// EXEC 


PGM= CALAHAN 




//FT06F001 


DD SYS<2>UT= A 




/ / FT05F001 


DD * 




Batched 


CALAHAN coded input (See Appendix E) 




/* 







Figure 11. Control Cards to Execute CALAHAN. 



The extensions and additions made to the CALAHAN program 

are: 

(1) An automatic printer-plot of the requested output for 
both frequency response, calculations and time response 
calculations . 

(2) A double-precision root-finding algorithm which 
usually works better than the original "Method of 
Mueller"* algorithm has been added. The roots are 
found by both methods, two polynomials are then con- 
structed from the two sets of roots, the magnitude of 
the coefficients are compared to the original. The set 
of coefficients with the least error determines the 
selected set of roots. 
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Where 



And 



V. 

Z( S ) = F(s) = 

in 



T(s) = F(s) 



V 

out 



V. 



in 



Considering now the case of non-ideal inductors and capacitors 

we can define a new function or M s M in which the impedances/ admittances 

become (r + Ls) / (g + Cs) where r and g are the dissipative 

O J_; O 

elements of the inductors and capacitors, respectively. Now defining 



S 






We have 



2,(5): L(s t sj and Y c (s) : C (s + S c ) 



At this point we could depart from the more rigoreous approach 
and make the assumption that ^ ^ ; thereby, making life 

very simple by drawing the conclusion that our function with non- ideal 
elements is formed by the substitution of (s + g ) for s in the original 
function. In order to test the validity of such an assumption the 
dissipation of capacitors and coils were measured using a General 
Radio Type 1650 impedance bridge. These components were chosen 
at random with the only criteria being that a ’’usage range” be obtained. 
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The results of these measurements are tabulated along with calcula- 



tions 


of % L and 


$ in table 
c 


2. As is 


apparent from the 


table the 




Type 


Nomial 


Measured D( 


) Cal. 


Comments 


Cap 


Electro 
. lytic 


4 ufd 


4 ufd 


270 


1.7 


Ignored 
see text 




Oil 


2 ufd 


1.93 


3.9 


0.0245 






Paper 


.0001 


.00093 


. 1.8 


0.0113 


Minimum 




Mica 


.0047 


.00446 


4.7 


0.0295 






Mica 


200 pf 


182 pf 


4.0 


0.0251 






Cermic 


250 pf 


258 pf 


1.0 


0.0628 


Maximum 


Coils 


2- r-f 


700 uh 


710 uh 


280 


1.76 






1- r-f 


200 uh 


206 uh 


265 


1.66 


Minimum 




RF - Choke 




58 uh 


105 


6 . 60 






Air - Dux* 


— 


19 • 5uh 


. 36 


2.26 






Air - Dux * 




742 uh 


1.2 


7.55 


Maximum 




Locally 

wound 




25 . 2uh 


.52 


3.25 






* Air -Dux is a trade name of 
transmitter- to-antenna lead 


a type of 
coupling 


coil used for 



Table 2. Coil and Capacitor Measurements 



numerical value of S L or did not vary excessively over the “usage 

range” except for the single 4 ufd electrolytic capacitor. As electrolytic 
capacitors have relatively little use where the pole and zero positions 
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are critical, this particular capacitor was subsequently ignored. 
However, the approximation that £ must certainly be held invalid 

or at least suspect as these parameters differ by factors of 100. 

By utilizing the digital computer, the analysis of a particular 
circuit can be reduced to the coding of a few cards. Utilization of 
the CALAHAN digital computer program in accordance with the 
instructions in Appendix E provides the desired quantities with a 
minimum of druggery. 

The two circuits used for this problem were: 

(1) The Pi-Section of figure 12 which is often used to couple 

r 

a transmitter to an antenna, and 

(2) The Chebyshev filter of figure 13. 




Figure 12. Pi-Section Coupler 




Figure 13. Chebyshev Filter 




A 
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r c 


V 

out 


V. 

in 




ohms 


T(s) = V. 

in 


Z(s) = I. 

in 


.01 


100 


-1.01, -.505 + j 1 . 321 


-.5812, -.2194 + j 1 . 3065 


.05 


20 


-1.05, -.525 + j 1 . 313 


- . 6265, - .2368 + jl . 3038 


. 10 


10 


-1. 10, -.550 + j 1 . 303 


-.6830, -.2585 + jl.3000 


.50 


2.0 


-1.50, -.750 + j 1 . 200 


-1. 132, -.4342 + jl.2565 



Table 3, Pi-Section Pole Movement (r = 0, vary ^ ) 

l-i 



Since it was suspected that the capacitors would be substantially 
"more ideal 11 than the inductors in their effects on the circuit the first 
computations were to vary over the "typical range" and determine 

the functions 1 pole positions of the Pi-Section. These are shown in • 
Table 3. Next similar computations for the "typical" range 

were made. These results are shown in Table 4. 





T(s) 


V 

out 


V. 

in 


- r 

°l-~ L L 


1.0 


-1.0, 


-1.0 + j 1.4142 


5.0 


-1.0, 


-1.5858, -4.414 


10.0 


-1.0, 


-1.228, -9.772 



Table 4. Pi-Section Poles vs. 



As can be seen from these calculations the effects of the typical 
inductor on the movement of the poles is some 100 times greater than 
the typical capacitor. This correlates nicely with the measured 
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The computations for varying £ c in the Chebyshev filter 
of figure 13 were carried out but the effects were again on the order 
of 2% (See table 5). 





Zeros 


Poles 


0.01 


-0.005 + j 1 .514 


-0.9819 






-0.2237 + j 1 - 140 






-1.217 + jl.425 


0.05 


-0.025 + j 1 . 5136 


-1.041 






-0.2446 + j 1.1378 






-1.2266 + j 1 . 4032 


0.10 


-0.05 + j 1 . 5 1 3 


- 1.116 






-0.2711 + j 1 . 1 356 






-1.2375 + j 1 - 3757 



Table 5. Chebyshev Filter (Poles and Zeros vs ). 




The transfer function of figure 14 was calculated and the poles 
determined for various values of RC1 and RC2. The particular values 
were 10, 100, and 1000 ohms which corresponded to a of 0 . 1 , 
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0.01, and 0.001, respectfully; thereby, bracketing the measured 

typical values. With each of these values RL was varied over a very 
-7 

wide range of 10 to 100 and pole positions calculated. These results 
are tabulated in table 6 and graphed in figure 15. 

As indicated by the graphs the pole positions certainly do not 
move straight in the negative real direction. In fact, there is quite 
a severe vertical movement. However, notice that the poles do come 
back down to level with the ideal at RL ~ 2 which is in the midrange 
of the measured ; thereby, in some sense justifying the n direct 
left shift 11 as a ball park design criteria. 

The effects of the capacitor admittance is to move the pole 
position slightly to the left and downward while the inductive resistance 
moved the pole up until RL L (numerically) and then the effect be- 
comes quite the opposite. The major point to be gleaned from this 
problem is that the usually neglected inductor resistance (in filter 
design) can be of such a magnitude to cloud the analysis (with ideal 
elements) in doubt. The CADLIB provides a tool to perform the 
mathematical druggery so that these parasitic elements can be included 
in the original design. 

ON COMPENSATIONS IN INTEGRATED CIRCUIT DESIGN 

Integrated circuits has become one of the most important areas 
in space and military electronics. As a consequence, ingegrated 
circuit design has become one of the major branches of network 
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RL 


RC1=RC2=10 


RC1=RC 2=100 


RC1=RC 2=1000 


io' 7 


-.5500 + j . 8352 


-.5050 + j . 8631 


-.5005 + j . 8657 


10" 6 


-.5500 + j . 8352 


-.5050 + j .8631 


-.5005 + j .8657 


10" 5 


-.5500 + j . 8352 


-.5050 + j . 8631 


-.5005 + j . 8657 


0.1 


- . 6000 + j . 8660 


-.5550 + j . 8905 


-.5505 + j . 8928 


0.2 


-.6500 + j .8930 


-.6050 + j .9143 


-.6005 + j .9163 


0.4 


-.7500 + j .9368 


-.7050 + j .9524 


-.7005 + j .9558 


0.5 


-.8000 + j .9539 


-.7550 + j .9669 


-.7505 + j .9681 


0.7 


-.9000 + j .9798 


-.8550 + j .9879 


-.8505 + j .9886 


0.9 


-1.000 + j .9950 


-.9550 + j .9985 


-.9505 + j .9987 


1.0 


-1.050 + j .9987 


-1.005 + jl.000 


-1.000 + jl.000 


1.1 


-1.100 + jl.000 


-1.055 + j .9990 


-1.050 + j .9988 


1.2 


-1.150 + j .9987 


-1.105 + j .9955 


-1.100 + j .9950 


1.3 


-1.200 + j .9950 


-1.155 + j .9895 


-1.150 + j .9888 


1.4 


-1.250 + j .9887 


-1.205 + j .9808 


-1.200 + j .9799 


1.5 


-1.300 + j .9798 


-1.255 + j .9695 


-1.250 + j .9684 


1.6 


-1.350 + j .9682 


-1.305 + j .9555 


-1.300 + j .9541 


1.7 


-1.400 + j . 9539 


-1.355 + j .9386 


-1.350 + j .9369 


1.8 


-1.450 + j .9368 


-1.405 + j .9187 


-1.400 + j .9167 


1.9 


-1.500 + j .9165 


-1.455 + j .8955 


-1.450 + j . 8933 


2.0 


-1.550 + j .8930 


-1.505 + j . 8689 


-1.500 + j . 8663 


2.2 


-1.650 + j . 8352 


-1.605 + j . 8037 


-1.600 + j . 8004 


2.4 


-1.750 + j .7599 


-1.705 + j .7190 


-1.700 + j .7146 


2.6 


-1.850 + j . 6614 


-1.805 + j . 6066 


-1.800 + j .6007 


2.8 


-1.950 + j .5268 


-1.905 + j .4461 


-1.900 + j .4369 


3.0 


-2.050 + j . 3123 


-2.005 + j .0999 


-2.000 + j .0316 


3.2 


-2.470; -1.830 


-2.555; -1.659 


-2.558; -1.643 


10.0 


-9.886; -1.214 


-9.887; -1.123 


-9.887; -1.114 


50.0 


-49.98; -1.120 


-49.98; -1.030 


-49.98; -1.021 


100. 


-99.99; -1.101 


-99.99; -1.020 


-99.99; -1.011 



Table 6. Non-Ideal Pi-Section Poles (Figure 14). 
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Figure 1^* 



Non- Id e til i’i- Sect ion ole 
(See Table 6) 



ft 
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synthesis. But solid state devices are extremely vulnerable to the 
effects of radiation, since the operation of most of these devices is 
determined by the density and gradients of the minority carriers. 

This minority carrier density can be drastically upset by radiation 
exciting electrons from the valence band to the conduction band of the 
material. When a transistor in a circuit such as figure 16 is exposed 
to a high energy radiation pulse the resultant collector current wave- 
form is similar to that given in figure 17. 

The integrated circuit transistor operating in the active region 
can be considered basically as one forward-biased and one reversed- 
biased diode connected back-to-back. When irradiated electron-hole 
pairs are generated throughout the semiconductor body including the 
substrate (e.g. , silicon). In the transition region near a reverse- 
biased p-n junction the generated minority carriers are swept across 
the junction by the electric field within a few nanoseconds. Carriers 
generated outside of the transition region will raise the minority 
carrier concentration gradient near the junction generating a diffusion 
current. These two components of current constitute the primary 
photocurrent in a p-n junction and; thereby, a transistor. The general 
shape of the primary photocurrent is very much dependent on the 
minority carrier lifetime as indicated in figure 18. In an npn transistor 
the generated collector region minority holes flow into the base region. 
This hole current is amplified by the current gain of the transistor 
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A/VW 




Figure 16. A dimple Transistor Circuit 



Dose 



Radiation Pulse 



time 




Figure 17. Transistor Photocurrents (induced) 
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Relative Primary 




Figure 18. Relative Shanes of the Diffusion 

Component of Primary Photocurrent. 
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circuit generating a secondary photocurrent such as that shown in 
figure 17. 

Radiation-induced photocurrents can cause errors and failures 
in integrated circuits. If the dosage is large enough the semiconductor 
material can suffer permanent displacements effects in addition to 
the ionization effects which are usually temporary [" 27^] . These 

effects as well as the detailed explanation of the theory and relations 
between the radiation dosage (electron, gamma, etc.) and the resultant 
photocurrent magnitude and shape we leave to those specializing in 
such matters. Of major concern here is the. fact that we have 
"parasitic" photocurrents whose effects must be at least diminished 
if not eliminated. One method is circumvention, which in its most 
common form amounts to removing power from critical circuits 
during the radiation pulse . However, in many applications (e.g. , 
counters, guidance mechanism's) continuous operation is mandatory. 




Figure 19* A Simple Logic Gate. 
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Dielectric (e.g. , silicon dioxide) isolation is a very effective 
means of removing many of the "evils" of radiation. The simple 
logic gate shown in figure 12 is used to explain the effects of dielectric 




Figure 20. Logic Gate with Parasitic Elements 

isolation. If this circuit is built on a silicon chip by a diffused process 
with p-n junction isolation (i. e. , isolation between the components 
and the substrate is obtained by reverse-biasing the p-n junction 
boundaries), the parasitic diodes and transistor shown in figure 20 
will generate substantial photocurrents making analysis extremely 
difficult. These photocurrents will vary widely and some will be 
quite large. Also the parasitic transistor forms an SCR and problems 
of latch-up are encountered. By use of thin film resistors and di- 
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electric isolation that is relative insensitive to radiation (wide for- 



bidden region) these parasitic elements and their leakage currents 
can be eliminated and the integrated circuit can be analysed very 
much the same as discrete circuits with, of course, the appropriate 
parameters. The resultant equivalent circuit is again figure 19. 



Figure 21. Photocurrent Generator and Compensation 

Diode (dotted) . 

The major problem is now the transistor photocurrent, i , 

PP 

as shown in figure 21. This photocurrent can be compensated for by 
adding a diode as shown (dotted). If thedode junction is made to match 
the collector-base junction, the photocurrent, i of the compensating 

diode will be the same as i so that no secondary photocurrent will 



PP 

result. The disturbance due to the transient radiation will then be 



0V CC = /0v 





o — KH-C4 




kept to a minimum. 



e 




Figure 22. Compensation Definitions. 



By defining the component compensation as 



1 N - (I - I ) 

% compensation = — — ( 100% ) 

N ' pp 



where the terms are defined in figure 22 it is possible to have some 
measure of the effectiveness of a compensation scheme at eliminating 
the effects of the secondary photocurrents. In the circuit of figure 21 
it is easily seen that at least some compensation will result if any 
compensating diode is used regardless of whether or not the collector- 
base junction is identical to the diode junction. Also since the photo- 
current of a p-n junction is directly proportional to the junction area 

f 25l the ratios of i and i , will be the same as the ratios of their 

PP PPd 

areas with the amount of compensation being proportional to this 
ratio . 
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In the circuit of figure 21 the amount of current compensation 

will be determined by how closely the junctions can be made to match. 

However, even with 100% compensation this configuration suffers 

from the voltage drop of i through the collector load which, in itself, 

PP 

may very well be intolerable. Also as pointed out by Grannemann, 
et. al. , £26^] the amount of compensation provided is very sensitive 
to a close match of the compensating diode and the transistor collector- 
base junction. This sensitivity is due primarily to the amplification 
of any mismatch by the gain of the stage. Utilizing the CIRCUS 
computer program £ 12 & 13^} an analysis of figure 21 was made 
using the parameters of the 1N903 and 2N404 as M typical n while the 
compensation diodes were special formulations to match (or mismatch) 
the collector junction of the 2N404. This analysis showed that even 
a primary photocurrent of 0 . 3 milliamperes peak was intolerable in 
that the output which was 4.23 volts with zero input dropped to about 
zero volts and stayed there for about four microseconds when subjected 
to a 10 10 rad per second 200 nanosecond radiation pulse such as shown 
in figure 17 . 

By attaching an additional identical compensation diode across 
the load resistor as shown in figure 23 we were able to achieve sub- 
stantially better compensation in that V remained more nearly 

out 

constant. The sensitivity of the output to a mismatch of this load 
diode is much less than the base as there is no gain following it. 

But as noted above imperfect matching of the base diode to the transistor 



53 




Figure 23. Base and Load Diode Compensation 

is magnified. Next over-compensating with the load diode was tried. 
This led to an ’’increasing" output voltage (with zero input) during a 
pulse. This factor was immediately seized on as a means to decrease 
the sensitivity of the circuit to the amplified mismatch of the base 
diode. The amplified mismatch lowers V ^ while the over compensated 
load diode raised . This over- compensation can be effectively 

accomplished by increasing junction area of the load diode. By 
utilizing a factor of 2 and our values a compensation of 80% was 
attainable for a base diode tolerance of 20%. 

Next the circuit of figure 24 was attacked. Initially only base 
diodes on T1 and T2 were utilized and theoretical 100% compensation 
was attainable for each transistor if the diodes and transistors were 
exact matches. However, sensitivity to mismatch worsened as T1 



54 





55 



mismatch is amplified by T2 when T2 is turned on by its own mis- 
match. As a tolerance of less than 1-5% is not feasible in integrated 
circuits with the present technology, this approach was abandoned. 

As T2 forms part of the emitter circuit of T1 attempts to com- 
pensate solely T2 met with very little success. By attaching base 
diodes to both T1 and T2 along with load diodes across the 1.75 K and 
the 6.0 Kohm resistors as shown in figure 25 adequate compensation 
was effected. As mismatch at T1 would cause a secondary photo- 
current into the base of T2 the base diode of T2 should be slightly 
greater in junction area than T1 base diode to provide maximum area 
matching tolerance. With the factor of two between base and load 
diode junction areas a compensation of 80% was attainable with areas 
tolerance of 20%. 

A theoretical development of an effective and economical 
means of radiation compensation in integrated circuits (also discrete 
but not so economical) has been demonstrated. An iterative transient 
analysis is required, thereby, virtually demanding the use of a 
computerized analysis. The costs associated with constructing, 
testing, and evaluation actual devices without the use of such an 
analysis would be prohibitive while with the analysis at least a reason- 
able first cut can be made. 



56 



BIBLIOGRAPHY 



1. Bashkow, P. R., "The A-matrix, New Network Description" , 
IRE Trans, on Circuit Theory, vol. CT-4, pp . 117-120, 

Sept, 1957 

2. Braddock, J. V., F. N. Hofer, and T. D. Guye, "NET-1R, 
The Radiation Effects Version of the NET-1 Circuit Analysis 
Program", Braddock, Dunn and McDonald, Inc., El Paso, 
Texas, Contract No. DA-49- 186- AMC- 365(X) 

3. Calahan, D. A., "Linear Network Analysis and Realization 

Digital Computer Programs: An Instruction Manual", Univ. 

of 111., Bulletin 572, vol. 62, No. 58, Feb, 1958 

4. Carpenter, R. M. and W. W. Happ, "Algorithms for the 
Dichotomous Representation of Systems", IEEE Trans, on 
Systems Sc Cybernetics, June, 1967, vol. SSC 3, No. 1 

5. Dawson, Darrow F. , Franklin F. Kuo, and Waldo G. Mag- 
nusom, "Computer-Aided Design of Electronic Circuits A 
User’s Viewpoint", Proc. IEEE, vol. 55, Nov, 1967, 

pp. 1947-1954 

6. Dervisoglu, A. , "State Models of Active RLC Networks", 
Coordinated Science Lab. , Univ. of 111. , Urbana, 111. , 

Rept. R-237, Dec, 1964 

7. Happ, W. W. and R. M. Carpenter, "NASAP: Present Cap- 
abilities of a Maintained Program", NASA Rept. ERC/ CQ 
67-615, NASA Electronic Research Center, Cambridge, Mass. 

8. Kuh, E. S. and R. A. Rohrer, The State Variable Approach 
to Network Analysis" , Proc. IEEE, vol. 53, July, 1965, pp. 
672-686 

9. Kuo, F. F., "Network Analysis by Digital Computer", Proc. , 
IEEE* vol. 54, June, 1966, pp^ 820-829 

10. Malmberg, A. F. and F . L. Cornwell, "NET-1 Network 
Analysis Program", Los Alamos Scientific Laboratory of 

the University of California, TID-4500, LA-3119 31st Edition 



57 



11. Mathers, H. W. , S. R. Sedore, and J. R. Sents, "Automated 
Digital Computer Program for Determining Responses of 
Electronic Circuits to Transient Nuclear Radiation ( SCEPTRE)", 
vol. I & II, IBM Space Guidance Center, Owego, N.Y., IBM 
File 66-928-61 I & II, AFWL-TR-66-126 vol. I & II, Feb. 1967 

12. Milliman, L. D., W. A. Massena, and R. H. Dickhaut, 

"CIRCUS - A Digital Computer Program for Transient 
Analysis of Electronic Circuits - User's Guide", The Boeing 
Company, Seattle, Was. , Harry Diamond Labs. Rept. 

AD-346-1, Jan, 1967 

13. Milliman, L. D., W. A. Massena, R. H. Dickhaut, and 
A. C. Mong, "CIRCUS - A Digital Computer Program for 
Transient Analysis of Electronic Circuits - Program Manual", 
The Boeing Company, Seattle, Wash. , Harry Diamond Labs. 
Rept. AD-346-2, Jan, 1967 

14. Pottle, C., "Comprehensive Active Network Analysis by 
Digital Computer - A State- Space Approach", Proc. 3rd 
Allerton Conf.. on Circuit and System Theory, Oct, 1965, 
pp. 659-668 

15. Seshu, S. and M. B. Reed, Linear Graphs and Electrical 
Networks, Addison- Wesley , Reading, Mass., 1961 

16. "1620 Electrotic Circuit Analysis Program (ECAP) User's 
Manual", IBM Application Program File H20-0170-1, 1965 

17. Sketoe, J. G. and Shu-Gar Chan, "ECAP/ 360 - A User's 
Guide, A Preliminary Report", Naval Postgraduate School, 
Monterey, Cal., Rept. NPS-52CD8041B , April, 1968 

18. Ludwig, A. , "Computer Programs for Antenna Feed System 

Design and Analysis - Volume I: Programs and Sample 

Cases", Jet Propulsion Lab, , Pasadena, Cal. , Tech Rept. 

No. 32-979, April 15, 1967 

19. L^udwig, A. , "Computer Programs for Antenna Feed System 

Design and Analysis - Volume II: Program Listings", Jet 

Propulsion Lab. , Pasadena, Cal., Tech. Rept. No. 32-979, 
April 15 , 1967 

20. "ECAP/ 360 - E: Electronic Circuit Analysis Program for 

the IBM System/ 360 Data Processing System", IBM S/ 360 
General Program Library 360D - 16.4.001, Aug 8, 1966 



58 



21. Decker, F „ R., W, C.Magnuson, Jr., A, W. Fulton, and 
R. H. Dickhaut, ’’Computer-Aided Design Review", Electro- 
Technology, Jan, 1968, pp. 74-76 

22. Bryant, P. R., "The Explicit Form of Bashkow's A Matrix 11 , 
IRE Transactions on Circuit Theory, vol. CT-9, Sept, 1962, 
pp . 303-306 

23. Wilson, R. L„, and W c A. Massena, "An Extension of Bryant- 
Bashkow A Matrix" , IEEE Transaction on Circuit Theory, 
CT-12, Mar, 1965, pp. 120-122 

24. Darlington, S. , Synthesis of Reactance 4-Poles, Journal of 
Mathematics and Physics, Vol XVIII, No. 4, Sept, 1939 

25. "Transient Radiation Effects on Electronics Handbook", 

DASA 1420, August, 1967 

26. Grannemann, W. W., et. al . , "Transient Radiation Effects 
on Microelectronics and Solid-State Devices", Technical 
Report No. AF WL- TR- 66- 1 29 , Jan, 1967 



59 



APPENDIX A 



CLASS "C" AMPLIFIER 

The following program selects appropriate parameter values 
to optimize a class "C M circuit. As rather complete comments are 
contained in the listing a detailed explanation is not included at this 
point, but two points should be brought out. The optimizing function 
is calculated in subroutine OPT and, therefore, it is easily modified 
to select or weigh some other parameters. Similiarly KCURR 
describes a tube model that can be quite easily changed to some 
other model (e.g. , a transistor model). 



60 





X X 


2 2 












<X <1 


«— *— 


2 2 










5 X 




k<— ' 










k— • k— < 


CC 


XX 










a cl 


a a 


<<] 














K h- 










• • 


• • 


UJUI 










l-UJ 


LL) 












o_j 


JC 


• • 




1- 






• • 


• • 


L— LLJ 




2 LL 


2 






hl- 


~<S 




a' o 


C* 




7Z 


xx 


• • 




O 2 


►— i 




i— *i— • 


CO 






UJ 


1- 




a or 


a ct 


<<r 




CL C 


< 








i — f— 




UJ CL 


a 




CL Ol 


a cl 


LL'U j 




a ll 


UJ 


e 










> 


K 


: z ’ 










1- —2 


k— 4 




LL LL 


LL UL 


U LL 




c cc 




<LU 


*— 1 1— i 


k— ) k*— « 


►-« *— < 




zee* o 


X 


cp 










• • G 


UJ 


LL 










2 C? O UJ CL 


c 


c: 


X 


2 


2 




CUJLUQCC 


LL 


2 


< 


»— » 


•— 




»-«cl cl a 






X 




X 




1 - oo 


2 


cf 


*— « 


C 


< 




OUOUOUOC 


C 


k— < »— » 


CL 


CL 


h- 




Coot/' LULU 




H* ♦— 






UJ 




CL LUlLOX^a 


<CL 










ULXXOC2X 


clo 






— * 


CL 


^OC J<IU 


LUS 


— * 


2 


2 


UJ CL 


JU-JU 


Ol 


2 


N»-» 


k— * 


XUJ 


>_j_j<ccr o 


OX 




5* 


X 


OJ 


o< <•— 




1- 


a 


O 


<1 


CL C 


ZHHhZhJ 


OJM 


3C 


CL 


K 


CL 


LUL-K— OLL-J 


x;* 


CL 


1 


LU 


1- 




1 1 k» i 


I- 2“ 


1 


K- 




Xi- 


UZZ-I 


LUO 


X 




1 


ao 


ll J k—« >— «t >— « w <J 




U’O — 


< 


O 




L-CL 


UlL^^ CL LL K- i 


mzh- 




cl 


<1 


X 2 LL' 


2LL UJUU 




•— <TO 


— o 


:JO 


1-0 


c*— CUJ 


C LULU LUO h- 


Cl 


S’ CZ 


a • 


a • 


uu • 


<o 


»- oo<~ 


O 1 


•-clo 


w 0 


— o 


— o 


UILUJ—C 




0-1 



hCU. 

cuuj 

<x 

O l- 



d 2 

LU CL »— 
2Uj*J 

~U_X 



II M 
< 



X 



h n 

CD 



II M 
O 



UJ 

* 

O 

O 

>■ 



-jjjh o_jhi--jll.-h 2»-<2 
aoa:-a-j>- <cjjoo *— llu i 
<<r>oo— a <oo> oclll-* 
L-l— >2Ct<K>> CLI-OlLUO 
aaiu tuooa lulu ►— 
ujuj|-OXr<\ lu>>i-ccoo<ujll 

OLX^O^KO J J<iK L~ LL 

uu - icluj xoclcl _jX«jc<lu 
< x< aorL xcl < o_ a a 2 x cc o 

LLH-L- OX < LU O- h— 
SJXI XU_r*j 2! 

D3w3qC*-C O XXlu2_juj 



UJ 

CP 

< 

I—* 

a: 

< 

> 



zxxo 
•-<>> 
S X X X 
O — <I<J 
CLOlXX 



2 

LU3- 

OCO 

LLLL 



COCO CO 



(X 

c 



oo 



cl cl cl cl 

CCOC 



UJUJULiUJ 

LLLLILU- 



<Lcl 

cc 



LULU 

LULL 



OO*— 1 1—. 






XXXXUJ3 X LU XX cc<u 




2 


oooo 


OO 


UjZ 






MMHMhOLU^hQMM 


ZD 




X 


<\J^vOCC 


CvJ^ 


CCL — 


o 




ZXXXZZ(£Z<mZX 


1— CL 


2 


X 


lilt 


1 1 


XX 


+ 




<1 <I< LL < X •— O CL •— < 


LL O UJ 


O 


— J 


rJ r— 4 r-4 rJ 




00c 


CO 




xx xxocoi-xacaux 


h-<q. 




O 


rovto 


CM 


oo 


+ 








♦— 


o 






r < 


<t 








CL 








xo 










►— « 








c: o 


ll 








CL 


f-l 






CL UJ 






2 




O 


V 






OU02 


o 




zxo-ui- — 2 




00 








CV)*- 






*~«<tl>> , 2 ►— X 




UJ 


o 


»H r— 4 rH H 


f\iC\, 


QCCLL 






XX XXLD -*X<CCOX 


£X<LL 


O 


CL 






OlJUJ 




o 


O •— < < O CP O »— CP O CC 


I-I-IL 




<r 






oo 




2 


aaXXLLLLOUJ>>>2 


2 LULL! 


h- 


O 






00 




LU 






X 








M 




o 






a 








X 




UJ 






2 








1— 




-J 






•— • 









(JOUUUOUUUOOOUOOUOUOUUUUUOOUUOUUUUUOUUUUUUUUUO 



61 



1- 


O 

or 


X 


►— 


l-C 


O 


X 


Ch 


“T - 


0 2 


O 


Zhw 


<T 


2d 


Hi 


on — d 


OO 


LL! oC 


2 or 


•— dl/' 


?o 


-J UJ 


a ll 


d<~ 


ThJ 


LUoO 


m<q. 


JQ 


ex 


Cl d 


c — • 


*-.<1 


LL 


H O 


d>0 


O 


O — d 


Xixi 


h-LU 


2 LU 


X<hJ 


d 


2d 1 


*-x 


<LL'2 


<r»- 



-Jl-C X 



U — 

oo 
> • 

o*-<Ln 

> 2:0 

XCOO 

<>*- 

►cr.— * 
crcr o 

>>«-• 
X ►f-H 



CSC 

I- 

2 



2 

l-X 


2 


CDi— • 2 


h- LL 

c? 

c 


<T<t — 
X 1-00 
►LUO 


x< 


CDOX 


Ol 


<*n on 


_J ► ► 


o»- 


CCOCD 


1 - 


h- 


►2~ 


O LU 


>>>22 


LULU 

Oco 


x*— • 0 

< d «— » 



o 

X_i 
< 



T3C - 1 
O.Q.— 

—■4 *>cC 



u 

> - 

XX 

<J 1—4 

5 2 
<? 

►— ► 

►uu2 

or 

> 

x 

< 

2 

►2 — 
h-»— ' O 
OlH 
XOcDm 

<io>cm 



o~ 

>2 

x< 



2 UJ2 



CC ►Z 
>h-~ 
XXX 
<CCD 

xo> 







O*-. 








3 ► 












UJK 


LUX ► 






H-4 sC 












17 


>0~ 






CL • 








CO 


OCQ 


O LU 


<dO 






O 


O 










1— x 


CC JpW 






CM 


O 




— 


dd 


ddd 


<IO 


♦ ' d • 






► O 






X 


OC 


COO 


CD LU 


*■ •> r 






2 ► ►CO 






< ► ► 




1—4 •— 4 


OO 


xzo 






— ZZO ► 


II 




X 3:0 


»■ •> 


— 


<1 


<— — 






X COO'v 






1-4 COO 


LULU 


LL LULU 


O 


23 ► 


— 




cll> ► 


2 




clll> 


— •» 




22 


o<— 


^2 




Q- O 








LL LL 


LL LL LL 


»— k-4 


~i-o 


1 — 


O 


• 


LL 


O 








O 


►LU--« 


1 X 




O 


O 


O 








GO — 


to ► *-* 


1 CD 


LU 


— CM 




r-l 








LU> 


ox — 


1 > 


X 


0-0 


LU 




— — 



00 00000 

vooc (Mvfsor^cc 

Mill 



I I 

r~W *— 4 

vTO 



osJnT vor- 



1-0 
<oc 
or cl 
LU 



ar<co x ► 
-r> >t x 
► * < 
OOCL — (MX 



Oh- st 

02 ► - 

sOUJ CD 

I OCP ► 



< 

> 



II 21 - 

* 

2CXLUCD 
COO CD 



r\j(M mmmmm 



CL> 


0 ►o mo 


* 


< 


QLL>0 




XI 


OCO 


OZh ^> 




I- 


2 


• 


1- 






mmH <— ■ - 


(T1HH 


< 


UJ 


O 


-J 


0 — 


XOJ 


►X— *-CM 


ao 


c 


► 


CM 




r-4 


<0 


oooomuj 


in • + 




0 


O 


< 


• O 


a: <r 


— -d >— <—o 


«-tO 


►— 


0 


M* 


LL 


LU CM 


OX 


►MVHlilZ 


>r 0 


0 


<M 




LU 


— J ► 


O 


CDO< — 2UJ 


»— » II O 


CL 


► 




C 


•sO 


d LU 


mVSOOwjo 


• X 


2 


in 






<w 


cl or 


2 20— J<I 


niL 


»— < 


w 


1— 


1— 


2 


< 


OOCC > II 


LUH 






< 


LU 


LL 


GO 


-JXX — -J — 


H a: 


c 


Q 


X 


00 


— K 


M 


<XX OQ 


xo 


< 


< 


OL 




k— 4 


X 


LUCO UJGQmmQ 


LU 


LU 


O 




LL QL 


h- 


dOO dLUXQ-l— X 


OC 


OC 


LL 








•— ♦ 






•—♦CM 




i 






H 






O 







r-xMm 



o 

CM 



OOOOOOOO UU OOOOO 



OOO OOO O 



62 



201 FORMAT (1H1, 40X, •***♦* INPUT DESIGN CRITERIA *',I3, 

0 • ******•///, 

1 9X » • POMI N = •♦1PEII.4,* WATTS • , AX t 

1 ' P I MAX = ' » E 1 1 .4 , * WATTS' ,4X, 'MAXVR = ',E11.4,' VOLTS', 4X 







C 


- 










► ► ► 




►-* 


00 










xvxx 




K 


0 














0 


u -_J 










► ► ► 




Z 


Z V.C 










• — — — 


•• 




*— 1 •» — ^ 










CtOCO 


*X 


LL 


ox- *x 










< a k* 


X. 




CO — — 










to cx O O 


X 


O 


► -0 










Q- <C C 


x^ 


Z 


0- o^<i 










OLL >> 


► 


*— i 


O C *fM 












• 


rxj 


00 >*— n 












* 




— r —4 — 










— — — — 


# 


•— 


GH - LL X 










#• — — 


* 


1 — 


0 < * -*x 










sf vf ^ 


* 


Q 


>^vf* * 










• • • 




c 


• >* 










«— < r— < — *■-< 






► — * — ^ • 










— *— 4 ^ H 




LL 


Z. — • — 4 II • 










U IL *LU 




0 


*■—’■' ^ LU •— 4 










— »•*— < *. 






? • CL LL' 










» — ^4 • 




UJ 


CC<M— <^CL 










UJ 


co 


n 












l| II CL || 


U 1 


—i 


LL* 3 C - 










*—• 


n 


< 


►Q- >* 










I- -z 


O 


> 


G *— Ml - 










»— * t— ' 


<r 




n - * n 










mo r - 


> 


_ j 


C- cnx 










on 11 cr 




<1 


CL' CC vt 0 






0 




a 0 :> 


LL 




15 11 > *0 






0 




V — — 


> 


h- 


CL - - - 






> 




x ► * 


t— 1 


►— • 


^ ► * 










<? XODX-* 


»- 


Z 


**— xoox 






0 


— 




c 


b—> 


CTCLCCI— vf 






z 


a: 


X -> *v 


or 




CO Z 5 h- * 




OOO 


c 


cc 




UJ 


LLJ 


>CL -< 3 - 




r\jr\!LP 




> 


OO co — 


h 


1 — 


• • ZT CO 


fx 




— 


c 


h- O* 




n 


• L- 


a 


one: 






OLO -O 




a 


ULXZ» -J 


LU 


khh 


*— r 


# 


oaxci 






LL LP LU *0 


2 




51 


z 


>O Cr >rv 




0 


LU *Osf> 




OOO 


a 


0 ~ 


• w- 




0 




Z 


coo 


> 


• T 


\ O 






*1— .oc: »— j- 


LU 






CV CD 


— * X * sC’ 




a 


2 ^ *LL' *— 4 * 


UJ 


«— 


— 


a: — Z > Z 






z 


►— « * LL LL 1 


cr 


z*— . 


cc 


Z CC -1 O -* 


t <r*x 


S' 


<1 


or • * • 






CP 


— > XZ # 2 


• • • • X 


UJ# 




^0* - — < 


< 




> 


57 O ac 0 cc 


r -4 *—i fH #— 4 X 


z# 


OL 


CLZ* * r-H 




Ooc>— 




CL O CC> • > 


r — 1 » 


»— * — 


LL! 


CXJ II UJ 


cc 




oc 


cr >0 # (ED (V D 


U-U II LL LO 


0 ► 


y~ 


<X *Q_ • * 


uu 


CL < 


0 


CD >CC CD> 


— — Q — « 


x —* 


z 


h- X LUf^O - 


K 


• h- 


a 


># CC O CD + 1 + 


*• ► 0 • * 


tr\^v 


— 


uj I— an 


fc — * 


a *uj 




> • >* 


• — 


— - 


0 


—00a 0 n 


cc 


)K~ 


00 


#LH #0 CM O Z Z Z 


II II - II 


rv» -vf 


0 


— # n. oor 


0 


• -J • 


0 


O — k M — » +~4 »— • 




0\<5 




0 O- -SO 




• UJ 


<x 


in •in+ac 1 o^ozocz 02. 


0 II 


(\l*X(Xj 


Z 


C? *C\J ► *CL 0 


LU 


h- O 


►— » 


000 0 0 ocooccaa: 




— x*. c< \ 


0 


w-o **x - - > 


UJ 


OX • 


H- 


• •a.w.ccwoL 1 >w>w> 


X ^ 2; O II 


"^3 **' —r*~ 1 


*—* 


hOCX ►- 


CL 


C <1 


a' 


O l| OCOh— CCH cLlL ►— H" 


<r x— 0 


— X + 


h- 


Q_r-»w— || X * 


X 


oc: 2 :< 


< 


>cl >cl><3 11 cl 11 a 11 


**x^> 




<r 


O h* cnx 


a 


JC»-I— 


a 


11 z 11 c 0 w 0 a 


- - < 1 - Z 


UJ< 


or 


0 11 o<a— 




CL CLUJ 




— II II II z z 2. 


h- • 


ax 11 


IL 


0 a? a 


0 


^ • — 


UJ 


CL 2 O O O II — O— iO—* 


JJ 


•— 'Ql 


h- 


II -J LL *— < QC LL! 


0 




h- 


aCCCOCQOCDOCC UUI 


» 


an 2 


»— 


<u csca- 


< 


a ll a 




>>>cc<cc<cn- 4 cu<cr.<rcr 




3U. Z 




^ouJ3:a 




>-— « »— « ► — < 


CL 


QCiC>U>U>D>U>U> 


<xj n >t m vC 




K- 


f-iruro 


CO 




X" 






CM 


LU 


ro 


< 


r- 


O 


O 




0 


OO 


0 


X 


•—4 


O 


CM 




rsj 




(NJ 











uu ouo ooo ooo 



63 



Uj 



5 Z 
□C O 
> 

+ X 

cc o 

> 7 

X 3 
< LL 
X 

— c 

Z 

# •— 
fsj 

to >— 

• X 

o l — ' 

X 

ii a 
o 



o 

Q 



O 

o 

> 



X 

ac 

> 



3 

o 

a 



tr 

a 

O'. 

► 

SO 



z 








GO 














w 










CO 


LL 








D 






X 








o> 


o 












ILf 


cc — 








> 










*> 




X 


> o 








X 


UJ 








cr. 




t— 1 


c o 








< — 


3 








tr 




CL 


> 








X CO 


-J 








> 


CL 


3: 


* o 








l/) II > 


< 










LU 












v* XX 


> 








•* 


X 




o * 




— » 




CM <J <1 










LL 


on 




• 




r . 


O' 


C2 X 5 


Z. 








LL 




to 


rsj O 




C 


o 


O' 


UJ 




O- 




LU 


X 


o 


• 




* 


> 


* > • 


Z 










X 


• 


( \ J 




ro 




oo on ►— 






O0 


O 


► 


LL 


o 






C 


c 


C 


LU 




X 


sfr 


*— 


►— 








+ 


z 


i— i < r°t *— * • 


1— 




-J 


tr: 


►<— < 


LL 


• 


z o x 




<\j 


< 


CX CO 


3 


<a 


3 


Ch 


QC 




X X 


*- o 




C! 




QC > — 


CL 


X 


<0 


x J 


Z3T 


X 


— J <H 


X > o 




* 


•- 


*- > *- xz 


X 


UJ 


UJ 


3 


CL 


QC 


• 


CC Q O 




cv 


z 


c <T~ 


a 




QC 


OOO 




LU 


o 


> O -It O > 




Q 


*— I 


CC OX 5 


o 






OLU 


•* 


Z> 


— o 


OOO 




+ 




cc + o cc 






X 


QC 




LU 


OL^OJ O 


* > • ># 






cc 


> > •> 


a 


* 


< 


— % 






•XO 'sf 


O r\J O 


2 : 


Q 


:> 


OZQh 1 


z 




•— i 


OX 




o 


LOiUJ <1 •* 




<r 






— O CD 


< 




Ci 


or 




UJ 


\ea — o 


-J + CC 1 — 14- — 1 


LLi 


r-4 


<► 


+ x ♦ • cc- 




O • 


UJ 


• <T 


o 


UJ 


uj<t z— io 


o a a a 


X 


c 


c r 


cc > 


s 


oo 


X 


c?cc 


r\j 


LL 


orv l xo 


| o— O — O 1 




««— 


cl 


CO O x rr 


X 


— o 


CL 


UiOJ 


•» 




^ <LZ LU •— < <r 


hohohuq: 


<1 


X 


> 


CO o<i o 


£ 


XX 


LU 


• CL 


n 


LU 


iC CX O'— 


3>Q->CL>C? 




ac 




> |l >5 • 




XCL 


1 — 


i 


w 


O 


II II <x 


— o o — 


UL* 


a 


x 


OM 


X 


+ e 


z 


QC U! 




< 


II II CL XX 


II II II 


X 


uQ 


2 


II Z II > — 


Z 


iC II 


•— 1 


XX 


LU 


a 


UJ LU — < 1 


II -J -J II 


< 




lU 


* — « w w 


UJ 






z 


X 




dJOUJO — xz 


O-JO-JO 


— J 


II 


X 


□0X0 


X 


II JIL 


o 


— C 


»— * 


UJ 


0<0<1 QC — 


f^u<iu<u^ 


3 




LU 


CCcCOlLLJL 


UJ 


<UL 


cc 


ULOC 


or 


> 


<CL<CLLLOU- 


o>o>o>o 


O 


oo 


* 


>>>^-^ 


CL 


iLO’JJ 


O 


•-*0 


3 


*— • 


CL Z CL Z LL *— 




-J 




O 




o 




O 


o 




o 






<r 




Z 




z 




LU 


LU 






LOO 




o 




*— 




t— I 




QC 


QC 






o<r 
























ro 


ooo ooo 


i OOO 


OOO 


1 oo ooo 


< 



ooo 



61r 



COMPUTE THE SECONDARY DESIGN RESULTS 





► 2_J ► ll 


on 


LI 


XyC< r II 


K 


2 


e <— 2 ooz 


1— 


► — < 


r- uuoock- — 


<r 


LU_J 


d ►—•2 ►►— <r 


3 


2 * 


► Xhl^U'X^C 




— H- 


OX or 03 




-JX 


i/" cl *— Q —* or 


X 


►c 


k- hi" - LLi 


X 


XCL 


•j oau X3 


► 


<3. ► 


C xcz II x X 






> m — o ►CL 


• 


<2 — 


sy <Tor >t 


CM 



— cc • & »-h 

- ► X^LU<0 3<M0 LL 











— 


X LU 


cc <?a at— c^l 


► 












<J ^ 


►cmT o oon luo 


X 










► 


2 — 


vfcr OX ►*-« 


cr 












OO 


• ^ Ok- _ ► 


► 










<r 


> •* 


(M>v ►rOXX || x 












<v 


LL 


r— "t x ►as: oc 


ii 










fO 


► UL 


ujoooxzca^ 












► 


LLUJ 


CL ^ m. — — LLJ ► 


a 












2 ► 


rH I ► ► OO 2 


LU 










► 


— CD 


►C X2 C • 


3 










* 


occ 


X > <1 1— CL CM 


C 










oo 


►or 


r- k* ►UJ 2 r- • 


o 










or 


X3 


► OvtCLCLOUJ 












UJ 


<tcl 


XO<TT2<t ► 


> 










k- 


2 


ll cr XMOhZx 


— J 










a- 


CL ► 


► r.ncm 


a 










2 


►— < LL) 


LLt ► ►h — H 


CL 










< 


►2 


cr- «xx xcvco ► 


X 










a 


X — 


<(MCCvO ► 


00 










CL 


<J 


K — 1 ► ►f-HXk* || 














2 ► 


o a <r ►in x 


LU 




. — ►— 






2 


CL CD 


cj ► • ►a 2 


k- 




o a 






cr 


>o 


> xojooooxk — 


<7 




oo 






*— • 


► 


oc^ahs Xor 


-J 




— . ^ xc • • 






CO 


►2 


V ►LL'S: O ►0 3 


a 




2 • OvJCV 






UJ 


LU — 


-J ►< O ^ X CL 


X 




— O O OJ 






— I 


2<l 


CL II ><MV 


r-( 




OO • • V Ss 








— o 


Ol CVfMk- 


rv 




— O (V (MV 


2 




-J 


o 


Olu II X c^ ►X 


► 




— — XX 


C5 




<T 


► ► 


oo cT ooisxc 


X 




I * * *CL*X < C 2 


— 






OLU 


<IU •‘Xvoa 


in 




a ->*“> <r 2 2 xx — 


O') 




— ' 


02 


UJkCvt ► ►'-< 2 


► 


— 02 


ox — — — x — 2 : xclclcc<oo or 


UJ 




LL 


> — 


k- o<r *Nt ►a 


X 




jj< o0002 2G<»OSa 3 — 


o 


LU 


• 


►O 


<Dk(\J • X 


X 


* >5T 


02 OO — ^>2 >CL3 a 2 




2 


► 


QO ► 


O > -J — t rv oo oo lTi 


•> 


O CG 


* cr **00200 o * * — a<. — 


V 


— 


X 


CC k- 


cl OiL'^ahn 




• 1 > 


“5> 2— — O ►O**- * k- x < 


QC 


-J 


o 


> X 


X>> ►ujik ► 


►<I 


f\J 


<r 2 mocncc xx wu o 


< 




CO 


a 


coo ►<<rx>t<M 


wx | 


W 4- V<1<UJ h- — 


o 


— 


► 


— 07 


CM CL -J 4U' 


• ro 


x<r 2 


^OCDwww<r 2 2 >XO * O 


2 


O 


—4 




► Cl <r l| ►CM ► 


— 2. CG2 


ooo * — * — cl cl cl 5 xaa«r oo 


o 


O 


X 


O CL 


x X 2 x Xr 




G OCO ► 


oua H 2 2 2 c<r~ > — ccjc •or — < 


o 


ro 


i-4 


CO ► 


ITOOOk II OOXVIL’V 




0>LX H — i. -Q.0f.O4O 


Li- 


► 


w 


►2 


► —2 »CD ► 


► 


il || — • Z>X oooooo OHM >300.0 


00 


o 




v£ — 


^C-OOIAJOK ►^ 


CD 


II II -> 


| II II O — — II II |l II > II CLH J 




W 


k~ 


— a: 


or o • nT <i 


O 


»— * 


-? ^ ->-> k-k-n ii <r 


k- 




< 


3 


k- or k* cr <x c\i *r\j 


► 


IXXin 


■ MwwXlLXZ II X Z>CO |l || 




UJ 


2 


UJCL 


<0 — DkHHfO 




Q, < < O 


ii -njHri/xx- CCcoh 2 u 


O 


k- 


CL 


k- 


s: xclo-Juj— i ► 


X 


JII 


— -o oo2 <2or2a:rt:a:X)<u.— 




k— « 


C 


— 


or cvlk* c ►LUX 




UJO£IC")^U<TCjU3:OCC'fl 'Jt^UlOhU. < CO 


k- 


oc 


LL 


or 


OmXj> X 




n>>o 


<ro>oo^— •►--•a^Q.aQ.QaiJLuoQ 


2 


3 




3 


UL 0< 








>— « 






«— t 


e-lfMCOstinNCk-* oo 


o o 


in o (mo 


QC 




O 




rH 




X vO 


vO o- r^ac 


CL 




a 




o 












CO 




CO 





uou o oo O 



65 



X 

o 



M 








•* 




X 






LO 


V* — > 




in 






O 


X. 




». A 






< 


•00 




s — 






a: 


LOO 




V V 






< 


UJ< 




M -w 






Ll 


•—•CL 




^ <5 






O 


3< 




< ► 






3 


ZU- 




cv >t 






O 


LU 0 




cn < 






►— » 


X3 




► <\j 








0 




x. ro 








MM 




Nv ► 






c 


XX 




► V 






3 


cr 




3 »“V 








M 




2X m 






3 


^X 




33 






3 


• ro 




org 






3 


OJ3 




a: mlo 






on 


H ► 




LU 3 : 






UJ 


3s* 




CL x 






ex 


CL • 




loO 








3f\J 




ir 






3 


►r-4 


X 


OCX 






3 


X LU 


< 


*ox 




— 


UJ 


O CL 


X 


>* co 




— CVj 


> 


r— 4 r-4 




• M 




* 


z 


M M 


LU 


fNJXsf 




— * 


0 


X 


3 


30C • 




3— 


0 


II co 


3 


3 ►rg 




CDH- 




M 


CL 


CL>*3 




ULZ 


3 


UJ 


X 


O *3 




LU 


z 


0 II 


O 


(\J •* 






<1 


z 


O 


r— 4 




X 




< LU 




LU II 




3 


LU 


KO 


• 


a 


3 


3# 


0 


OZ 


UJ 


H LU 


M 


O 


z 


3< 


0 


- -o 


3 


3*-* 


< 


Oh 


3 


z 


O 


a. 


3 


Zn- 


< 


II II < 


a: 


* ~ 


Ml 


MU 


3 


3 


MM 




O 


< 




>3L0 


O 


M# 


< 


— O CL 


c 


UO w 




CL 


CL 


K Z< 


0 


2200 


* 


O 


< 


3 0^0 


3 


UJ< UJ 


z 




O 


a z 




*—* 3 3 


< 


* 




00 00 


X 


OLO 


3 


O 


H 


1 0033Z 


3 


mm ro 




• O 


3 


0 * X*- 


O 


3L0< 


3 


eg • 


CL 


wo -02 


MM 


LL IllO 


3 


wvt 


3 


C0(NJCSJ3 




333 


CL 


w 


3 


* ►O *3 


O 


X 


►— 


x. 


C 


Orox O 


z 


(^hD 


3 


V. 




gJO *000 


<x 


hD* 


O 


O 


LU 


UJ vO-^UJO 


CO 


•*0.0 




•O 


co 


0 wwom 




XZLU 


LU 


»— 4 • 


3 


30 


0 


o»-mqc: 


z 


w 3 


3 


HOUJ«0 


UJ 


HII 


Ml 






-kxih 


3 


•'O' sf 


X 


If II 


LU 


11 —33 


LO 


HM 


3 




> 


uocoogo 


3 




UJ 


03 


O 


Om-^lL 0 


3 


O'O’-* 


3 




X 




a 




UJ 




UJ 


nOcvj 


3 




Q 




3 


coo 


or 



cn 

ooo uou uoo 



• 


vf 






0 


• M 




CC' 


3 


4*m 




MM 


3 


r-4 00 






3 


33 




• 


CL 


CL mm 




% 


X 


M^or 


»> 




O 


MX 




z 


mO 


• 3 


* 


0 


>* 


LOX 


X 


M- 


• 3 


MM 


r-4 


LO 


r\j-j 






3 


r-4 CO 


Xm 


O 


O 


3m 


3 M 


eg 




CLLO 


o>* 




3 


r-4 OO 


3 • 


X 


O 


mQ 




r- 


U 


• CL 


0-4 


M M 




U 


Z3 


l A* 


LO 


or 


<CL 


— < 3 


0 




CC r-4 


mO 


z 


11 


M 


X3 


0 


l—MM 


3* 


vC3 


0 


— OX 


33 


M> 


3 


(\j M< 


cco 


-*Z 


LO 


* JI 


<2 


* O 




* 0 


Z< 


TO 




— Z 3 


Mh 


— M 


— 


♦- <rx 


<0 


w 3 * 




— Z CQ3 


30 


300 


M 


1-3 


hC 


Zege* 


r\j 


30 O- 


<z 


1— 1 m3 


3 • 


03 3 - 


MM 


3X> 


3 s c 


O HX 


x 


MHZ 


3 3 


•w- LOO 


30 


3 mO 


r 


# 3 3 r-4 


3r z 


zvu 


M-* M 


M-.Z3 m 


Zmmmm 


mV 


h 


3 CL mm O 3 


MMXZ 


3 »*□ 


O II 


3^~33>v 


3<3 


M-h 


1 O 


c or •* 


mZK 


Z^ 


5L 3 


OC* •*• - 


3m 


M MQ 


— X 


3 **3 


M MX 


3>t 3 


O MM 


IfKCCXO 


JXK 


Z< 3 


wllJh 


3303 


QDOm 


MMfV MM 


3Z 


mQ 3 <3 


•*-13 


3rn < 


r mz 


30 • m 3 


3 M*> 


M MU.— 


—.h-O 


33 


2 V m^ 


3 X >vl- bT CL 


* * ZsO 


M-VXV 


2hl0VmO* 


mm r-O 


3 Min>v 


— ^Z — 


M 


003^Km >tH ► 


3(MO*st * X 


• • ^ ' 




mmm<J 


— O 


cgs*— r\jL03 3r\j3< 


— ^ K— OJ r-4 r-4 r-4 


— ' O^M\Om M(\j 


cs<mco m 


W O M 


MnT CO 


ova: • 


* *ov 



v* ^ mxlo mo inviuxoir.s 



V. ►KK'OHClXO *“M~h 

Q vOw fsjwwOHO'O^M •* || O'v 
• O^ CO ►LPlWr N* I 



— M 


• h- 


- V 3 


V 3 


•^3 


3 




r-»3< 


V.3<m 


V.03<- 


(\JX3<0 


II 


i-r 


3X 


KKX 


mm3 X33 




11 ^-3 


mm q: 


MM 3 


3mm CL 


OO 


'X 


30 


30 


030 


X300KZ 


033:3 


33 


033 


O3 3OL03 




r-4(^J r— 4 C\J -M(NJ 






0 


O 


r-4 


OO 


O 


r-4 


O 


O' 


O 


O 


OO 


O 


0 


0 




's* 




Olf\ 


O 


0 


0 








gr 


in 


in 


0 



66 
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FUNCTION SIMP ( A , N ) 

DIMENSION A ( HO) 

PERFORMS A SIMPSON RULE INTEGRATION OF VECTOR "A" 
NN = N+ 1 
A (NN+1) = 0.0 
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APPENDIX B 



FILEMAP AND FORMAT 
CIRCUS DISTRIBUTION TAPE 



MAP: 



File 1 CIRCUS OBJECT DECK: From FORTRAN H OPT 

0 COMPILER and release 9 assembler. 4000- 
5000 cards. May be loaded directly from tape. 

File 2 CIRCUS MODEL LIBRARY UPDATE PROGRAM 

OBJECT DECK (complete) approximately 1000 
cards . 

File 3 CIRCUS FORTRAN SOURCE DECK. Approximately 

6000 cards. 



File 4-12 CIRCUS ASSEMBLER LANGUAGE ROUTINES: 
Approximately 500 cards total. 

File 13 CIRCUS MODEL LIBRARY UPDATE PROGRAM 

FORTRAN SOURCE DECK, exclusive of those 
routines common to CIRCUS: See CIRCUS PRO- 

GRAM MANUAL (No. 346-2) less then 500 cards. 

File 14 CIRCUS SAMPLE DATA. Approximately 500 cards 



FORMAT: 

The 14 files fill approximately 1000 feet of tape. 

All files on all tapes are 80 character card images: 80 byte 

records unblocked. 

All tapes are odd parity, density 800 bits/inch. 

Note that the first two files only are in object form, all others 
are in EBCDIC. 



69 



NOTES ON CIRCUS/ 360: 



In order to CIRCUS/ 360 as supplied, the user must provide 
timing routines which are called by the program. These routines are 
a part of Lincoln Laboratory's system library and therefore are not 
a part of CIRCUS. Dummy routines may be written, but if different / 
names are used, or if similar routines exist in the user's system but 
with different names, the calling statements must be changed in 
CIRCUS. Calls appear in KLOCK and OUTPUT in several places in 
these routines. The routines are: 

TIMDAT (IYR, IDAY) 

TIMHR (HOURS) 

TIMHMS (IHR, IMIN, ISEC) 

where I implies fixed point and 

H implies floating point numbers 

These routines use the system clock and times as implied in their 
arguments . 

If the system on which CIRCUS/ 360 is to be used has a core 
capacity greater than 256 K bytes, the capacity of CIRCUS may be 
increased to handle circuits of greater than 70 equations by adding 
C0MM0N Al(50000) CIR25715 

and changing CIR257 50 to read 

LC0RE = 50000 CIR25750 

in SUBROUTINE MAPPER. 
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In order that the supplied version fit into a 256 K byte core, 
the system's load point must be B000 or less. 

CIRCUS/ 360 uses the system clock via the assembler language 
routine SETLMT to compare execution time against a specified limit. 
If no limit is specified, the program assumes a ten minute limit, 
after which CHARET is called and execution terminated. To adjust 
this limit, the first card in the data deck should read, beginning in 
column one , 

LIMITS, X 

where X is the desired limit in minutes, floating point (as in LIMITS, 
8.5 for 8 min. 30 sec. maximum execution). 

CIRCUS/ 360 now includes a printer plotter routine which is 
called by the PLOT statement in the input list (see user's manual). 
Plots are 51 by 120 characters; time interval is the smallest in the 
INTERVALS statement; if greater than 120 time points at this interval 
exist, the plot will be continued on succeeding pages. One plot vari- 
able appears per plot. 

KNOWN BUGS IN CIRCUS/ 360 

1. Occasionally a circuit time and associated outputs may be 

printed twice. This occurs at times when something is 
specified in the input list as changing. The first of these 
printed times represents a value slightly less than that 
printed by an amount determined by machine accuracy, and 
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rounded off when printed. The second value is the exact time 



as specified in the input list. 

2. When specifying a logical unit for a model library, save tape, 

etc. , a user is supposed to be able to specify the logical 
unit number (as in SAVE, 15) or use a letter/ number 
equivalent (as in SAVE, E) where the letter / number cor- 
respondence is specified in LINK2 (DATA MNSAVE/ 1HA, 

1HB , 1HC , JHD, 1HE , 1HF , 1 , 2 , 3 , 4, 1 5 , 1 6 / CIR20600). 

As it now stands, the number designation will not be 
recognized, while the letter designation works. If a number 
other than those in the standard letter/ number correspondence 
is desired, the aforementioned DATA statement may be 
changed to provide a more useful correspondence. 
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APPENDIX C 



THE CODING OF POTTLE 
LINEAR NETWORK ANALYSIS PROGRAM 

Given an arbitrary, consistent, linear, fixed electrical 
network (with zero initial conditions), the POTTLE program cal- 
culates: 

(a) . A prescribed but arbitrary transfer function t(s) or, 

more generally, an arbitrary n-port matrix t_(s). 

(b) . The poles and zeros of t(s) or t_(s)„ 

(c) . The frequency responses associated with t..(s) - 

prints and printer-plots magnitude, phase, and delay. 

(d) . A set of state equations for the network. 

In addition, any set of element values may be stepped up to 
six times and the above analysis repeated at each step. 

The network size limitations are: 

(a) . Maximum number of elements - 60 

(b) . Maximum number of LC elements - 32 

(c) . Maximum number of nodes - 27 

(d) . Number of controlled sources plus number of inputs 

must be less than or equal to 32. 

(e) . Number of controlled sources plus number of outputs 

plus number of resistors must be less than or equal to 64. 
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(f) . Three time (number of LC) plus two times (number of 

inputs) plus number of branches minus number of 
nodes must be less than or equal to 128. 

(g) . Maximum number of stepped element values - 88 



The input data format is: 




(a) . A title card - 


columns 1-80 will be printed on the 


output sheet. 





(b). Cards describing the branches of the network: 



Columns 


Des cription 


1 - 4 


Branch name (left justified) beginning 
with: 

V - voltage source 
I - current source 
R, L, or C - resistor, inductor, 
or capacitor, respectfully 


6 - 10 


Name of controlling branch (left justified) 
if col 1-4 contains a controlled 
source; otherwise leave blank 


11-12 


Positive node (right justified) 


14 - 15 


Negative node ( right justified) 


Column 


I ^ Column 



H-12 O O 14-15 

Number nodes consecutively from zero 
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Columns 


Description 


18 


If branch variable is to be considered 
an output quantity, use symbols V 
or I as appropriate. 


20 


Type of control (V or I) if controlled 
source; otherwise blank 


21 - 30 


Value of element or strength of con- 
trolled source; independent source 
leave blank. Must have a decimal 
point (FI 0.0 format) 


31 - 40 


Values of stepped elements or 

strengths of stepped controlled 


41 - 50 


sources; all stepped elements 
must be stepped the same number 
of times; stepped element val- 
ues may not be continued on 


71 - 80 


succeeding cards* Must have 
decimal point (F10.0 format) 



(c). A card denoting the types of responses desired* This 
card must be present* 



columns 11 ■ 


■ 12: 01 if frequency response is 

desired 




00 or blank otherwise 



(d). Two cards to describe frequency response desired; these 
are omitted if no frequency response is desired. 



Card 1: 




column 1: 


1 if linear frequency scale 

2 if log frequency scale 



column 3: 1 if magnitude calculated in db 

2 if magnitude in ohms 

3 if magnitude in mhos 
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column 



linear frequency scale 

total number of frequencies 

5 - 7: 

logrithmic frequency scale 

number of frequencies per 
decade 

number is right justified 
(13 format) 

Card 2 

columns 1 - 10: Lowest frequency (Hz) 
columns 11 - 20: Highest frequency (Hz) 

(e). If more than one network is to be analyzed, a set of these 
cards is required for each additional network. 

Additional comments: 

(a) . The method of analysis used is sensitive to scaling, 

producing overflows if improperly scaled. Values for 
L’s and C's should be within an order of magnitude of 
unity for large networks. 

(b) . Zero values for inductors and capacitors may be used 

to establish a dummy branch current or node-pair 
voltage; this is useful in controlling a controlled source 
or defining an output variable. 
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APPENDIX D 





THE POTTOL TRANSLATOR 


PURPOSE 





The primary purpose of the Translator is to convert an 
engineer- oriented free-format circuit description into the standard 
input format required by C. Pottle* s Linear Network Analysis 
Program. Prior exposure to the POTTLE program (See Appendix 
C) is assumed. 

FEATURES 

The major features of the Translator are: 



(1). 


Free-format of input data 


(2). 


Minimal changes to the original program 


(3). 


Simple and effective overlay (if needed) 


(4). 


Additional diagnostics 


(5). 


Exact reprint of input data 


(6). 


Continuation cards 


(7). 


Comment cards 


(8). 


Batch processing 


OVERLAY 





The particular form in which the Translator has been con- 
structed lends itself to a very natural and fast overlay structure. 

Figure D-l shows the recommended overlay structure. This particular 
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Figure D- 1 . Overlay Tree Structure. 



structure loads the Translator along with the root segment, processes 
the input data, loads the rest of the program, executes the program 
and stops without having to repetitiously load and reload the same 
overlay branch. 

POTTLE PROGRAM CHANGES 

Since the logical unit assignments used in the original program 
were not the same as used locally, it was judged worthwhile to change 
all logical unit references to variables and assign them at one central 
location. This is presently accomplished in the main program of the 
program. These assignments are transferred via labeled COMMON 
to the various subroutines. The present assignments and their 
correlation with the original assignments are shown in figure D-2. 
These assignments may be changed permanently by changing the 
appropriate cards in the main program. Additionally a temporary 
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Variable 


Unit 


Use 


Original 


INTAPE 


5 


User's input data 


i 


IOUTP 


6 


Copy of input and errors 


2 


STATP 


6 


State equations 


3 


TRAN TP 


6 


Transfer functions, poles/ zeros 


4 


FREQTP 


6 


Frequency response 


5 


NOTE: Unit 6 is the printer locally 



Figure D-2. I/O Unit Assignments. 



change of any or all of units may be invoked for a single run through 
the use of a TAPES card as shown in figure D-3. 



column 7 

X 


or later 




TAPES, 


STATP = 2, 


TRANTP = 3 



Figure D-3. Temporary Unit Assignments 



ADDITIONS 

The additional features added to the allowable input include: 

(1). Continuation Cards: Any number of continuations are 
allowed as long as the previous card's last character 
is a comma. There is a limit of 154 non-blank charactors 
after column 6 of the first card. 
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(2) , Debugging: A limited debugging facility is provided in 

that a card containing the word TRACE (starting in 
column 7 or later) will cause a diagnostic print as control 
enters and departs each subroutine. It also provides a 
print out of the NHEX variables in A1 format. 

(3) . Batch Processing: Batch processing is implemented 

by requiring that each set of circuit description data 
end with an EXECUTE card and the batch run or job 
end with an END card. 

(4) . Comment Cards: Except for reproduction in the print 

out, any card with an asterick in column one is ignored. 

(5) . Number Format: All numbers on the input data cards 

may be either F, E, or integer format with or without 
decimals. Recognizable number forms are: 

-XX. X E-YY XX . XXEYY 

XXXE+YY X.X + .XX 

If the decimal is missing it is assumed to follow the 
last X. 

(6) . A print out of the translated program is provided. 

INPUT DESCRIPTION 

The standard R, L, C, V, or I branch is coded by placing a 
unique name for the branch in columns 1-5 (note: This is the only 
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type of card allowed to have anything in the first five columns). If 
five characters are used in the name the fifth character will be 
ignored; however, the name need not start in column one. All other 
information is placed after column six (including continued cards). 

The node numbers are delimited by parens and separated by a comma. 
The order of the node numbers implies the positive direction of 
current flow. Spaces may be used freely to enhance readability. 
Typical examples are shown in figure D-4. In the case of stepped 



1 


2 t 3 ! 4| 3 


6 


7 | e I 9 I 10 


M | 12 | 13 | 14 1 13 | 16 | 17 | 10 | 19 | 20 
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1 \Z\0\ 1 1 


1 1 1 1 1 1 1 1 1 


V 


, , ■ 




,(A 


ll^lllpl 1111 


1 1 1 1 1 1 1 1 1 


c 
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) 


\C\ i7i)i l 1 


1 1 1 1 1 1 1 1 1 








> 







Figure D-4. Input Data Examples. 



elements the initial value (or only if not stepped) is placed first with 
sequential stepped parameter values separated by commas following. 
All stepped parameters must be stepped the same number of times. 

An output of the branch voltage or current is indicated by adding 
the expression OUTPUT = V or OUTPUT - I, respectfully, to the 
branch description. 

Descriptions of controlled sources have one additional restric- 
tion. The controlling branch name must appear BEFORE the type 
of control. This is illustrated in figure D-5. 
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Figure D-5. Description of Controlled Sources 
FREQUENCY RESPONSE 



The inclusion of a card of the form shown in figure D-6 will 
generate the proper commands to obtain the indicated frequency 
calculations. This card will set the variables via labelled common 
required in subroutine TPRINT. It is not necessary to include the 
words LOG and/or DB if these options are desired. 

COMMANDS 

The EXECUTE card causes an analysis to begin. It also 
generates the blank card required by the POTTLE program if no 
frequency card is included. 



FREQ = nl (n2) n3, LINEAR or LOG, OHMS or MHOS or DB 



nl-min 


freq (Hz 


Defaults 


n2- no . 


freq/ decade (log) 


FREQ 


none 


no. 


freq (linear 


freq scale 


LOG 


n3-max freq (Hz) 


units 


DB 



Figure D-6. The Frequency Card. 



82 



An END card terminates the particular run. Lack of an END 
card does not adversely affect the execution of the program (if all 
other data is correct). 

COMMENTS 

The Translator is intimately tied to the POTTLE program. 
Therefore a proper appreciation of the Translator requires at least 
a passing acquaintance with the use of the original POTTLE program. 
To this end Appendix C contains a description of the original program 
and the CADLIB contains the original program. 

In order that POTTLE changes could be kept to an absolute 
minimum the Translator operates on only one card (and its con- 
tinuations) at a time, i.e. , the order of input is the same as the 
original POTTLE program except that the label card is no longer used 
as unlimited comment cards are available for expositor comments. 
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APPENDIX E 



CALAHAN - A USER'S GUIDE 

CIRCUIT LIMITATIONS 

(1) Complexity 

(a) . 30 nodes 

(b) . 100 passive elements 

(c) . 20 active elements 

(2) Voltage controlled current source is the ONLY 
allowable active element. 

(3) Storage - approximately 90 K bytes 
DEFINITIONS 

(1) Network - the circuit to be analyzed when considered 
as a two-port black box (See figure E- 1 below). 



^ O 

^ Output 

— o 

Figure E- 1 . The Calahan Network. 

(2) Network function - A functional relationship between the 

input voltage (V ), the output voltage (V ), the input 
1 £ 

current (I ) , and/or the output current (I ) where these 
terms are positive as indicated in figure E- 1 . 



Input 
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(3) KEY 1 - A computer flag whose value indicates a particular 
network function (See Table E-l below). 



KEY 1 


Network Function 


Symbolic 


1 


Voltage Transfer Function 


Yl 

v, 


= 0 


2 


Open Circuit Driving 


v, I 






Point Impedance 


1 


1 J 2=0 


3 


Open Circuit Transfer 


V, 






Impedance 


- u. 

I, 




4 


Short Circuit Driving 


X, 






Point Impedance 


v, 




5 


Short Circuit Transfer 


I 2 






Admittance 


V, 


V 2 -0 


6 


Current Transfer Function 


I* 








X, 


V 2 =0 



Table E-l. Definitions of KEY 1. 



(4) KEY 2 - A computer flag whose value indicates a particular 
set of desired outputs (See table E-2). 



INPUTS 

(1) Network elements or network function (Voltage or current 
transfer function only. 

(2) Input and output node numbers 

(3) Type of network function desired as output 
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KEY 2 


Output(s) is (are) of the Network Function 


1 


Coefficients, poles, and zeros 


2 


Coefficients, poles, Zeros, and freq response 


3 


Coefficients, poles, Zeros, and time response 


4 


Coefficients, poles, zeros, freq response , and 
time response 


5 


Symbolic form of the network function 



Table E-2. Definitions of KEY 2. 

OUTPUTS (See Table E-2) 

(1) Coefficients of requested network function 

(2) Poles and zeros of the requested network function 

(3) Frequency response (both tables and printer plots of 
magnitude, phase and delay) 

(4) Time response to an arbitrary input (combination step 
in CADLIB version). Tables and printer plots are 
included, 

(5) Repeated outputs of 1-4 with some element value as a 
variable parameter 

(6) Designated network function in symbolic form 
INPUT CODING 

The initial step in coding for CALAHAN is to number all nodes 
sequentially from one. A deck of cards is then coded as indicated on 
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Figura E-2. CALAHAN Coding Flowchart 







the flowchart of figure E-2. Two points need clarification. The first 
card is a title card and is ignored except for titling the output. The 
second card is a dual purpose card. If the input is in network 
element form, card #2 contains a positive number equal to number 
of circuit descriptions between title cards. As each circuit or circuit 
variation will normally be different the usual value will be one. If 
the input is to be a transfer function the value used will be a negative 
KEY 2. 

DISCUSSION 

General Precautions 

(1) . Poles, zeros, and time response may be impossible to 

find for no predictable reason. 

(2) . Some bandpass networks may yield incorrect responses 

due to single precision limitations. 

(3) . The nodes must be numbered 1 through N, the number 

of nodes; otherwise, there are no constraints on the 
manner of numbering. 

(4) . The network must be connected (i.e. , a path must exist 

between any two nodes). Thus, for a transformer, the 
common connection between windings must be made, 
unless a path connecting the windings already exists 
through external circuitry. 
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(5) . The RLC elements must be listed before the controlled 

sources, although in any order. 

(6) . Two or more RLC elements across the same nodes must 

have nodes listed in the same order. 

(7) . Negative RLC elements are permitted. 

(8) . A zero value for a resistance or inductance is never 

permitted. 

(9) . If repeated analyses are made with some element value 

as parameter, a zero value is never permitted for a 
resistance or inductance, and is permitted only as a 
final value for a capacitance or controlled source. 

(10) . When a logrithmic frequency scale is used, recall that 

zero frequency is not a valid lowest frequency. 

(11) . When choosing a highest frequency, beware of possible 

overflow. If a^ is the coefficient of the highest power 
(n) of s in the denominator, then overflow may occur if 

2 n 7 g 

(2 Jt f a ) is greater than ( 10) 
n 

(12) . The time response cannot be calculated for a network 

function with a numerator of higher degree than the 
denominator . 

(13) . The topological program can be fooled by certain de- 

generate cases which do not warrant safeguarding. 
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